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THEORETICAL  STUDIES  ON  FREE  ELECTRON  LASERS 


The  subject  of  this  investigation  has  been  the  Free 
9- 

Electron  Laser'i"(FEL) ,  a  device  which  is  capable  of 
converting  the  kinetic  energy  of  a  relativistic  electron  beam 
into  coherent  electromagnetic  radiation.  The  investigation 
focused  on  two  main  topics: 

-A-r- '  The  FEL  operated  as  an  oscillator  and  an  amplifier 
using  variable  parameter  wigglers. 

"  The  FEL  oscillator  operated  in  conjunction  with  a 
storage  ring  using  gain-expanded3  and  phase  area 
displacement  wigglers.^  - - 

During  the  two-year  period  extending  from  August  1, 

1981  to  September  30,  1983,  the  following  tasks  were 
accompli shed : 

1.  An  eigenmode  analysis^  of  the  linearized  one¬ 
dimensional  FEL  equations  leading  to  a  determination  of  the 
linear  gain  characteristics  of  FEL  oscillators  with  variable 
parameter  wigglers.  Various  cases  studied  include  electron 
pulses,  long  and  short  compared  to  slippage  length;  constant 
and  variable  parameter  wigglers  ;  and  FELs  with  and  without 
optical  sideband  suppression. 


2.  The  development  of  a  one-dimensional  self- 
consistent  particle  and  electromagnetic  wave  code^  to 
simulate  the  operation  of  FELs  with  variable  parameter 


wigglers  and  with  frequency  filtering  of  the  electromagnetic 
(EM)  pulse. 

3.  A  study  of  the  time  evolution  of  an  EM  pulse  in  FEL 
oscillators  in  which  the  growth  of  the  pulse  is  followed  from 
low  noise  levels  all  the  way  to  steady  state,  with  particular 
attention  paid  to  the  effectiveness  of  electron  trapping  and 
to  the  suppression  of  sideband  modes  by  frequency 
discrimination . 

4.  A  study  of  the  effect  of  unstable  sideband  modes  on 
high-extraction  FEL  amplifiers,  with  particular  reference  to 
a  proposed  ATA  design.5 

5.  An  extension  of  the  analytic  linear  theory  of 
unstable  sideband  modes  to  include  not  only  the  amplitude 
perturbations  previously  considered, ^  but  also  the  phase 
perturbations  of  the  EM  pulse. 

6.  A  Hamiltonian  formulation  of  the  electron  equations 
of  motion  for  the  "thin  lens"  gain-expanded  FEL, 6  with 
derivation  of  the  Manley-Rowe  relation  and  generalized 
gain-spread  theorems. 

7.  A  study  of  the  phase  area  displacement  wiggler, 
with  estimation  of  the  theoretical  efficiency  possible  in 
storage  rings. 
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8.  A  preliminary  survey  of  the  role  of  two-dimensional 


effects  on  FEL  operation. 

Tasks  1,  2  and  3  for  ultra-short  electron  micropulses 
as  well  as  most  of  Task  6  were  completed  by  the  end  of  the 
first  contract  year  and  are  discussed  in  detail  in  the  1982 

4 

Annual  Report. 
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I.  FEL  OSCILLATOR 


In  an  FEL  oscillator,  a  variable  parameter  wiggler  is 
positioned  between  two  mirrors  and  an  EM  pulse  is  reflected 
repeatedly  through  the  wiggler.  Relativistic  electron 
micropulses  are  injected  at  periodic  intervals  so  that  on 
each  forward  pass  of  the  pulse  through  the  wiggler,  there  is 
overlap  of  the  EM  pulse  and  the  electron  micropulse.  The 
pulse  grows  in  amplitude  on  each  forward  pass  and  eventually 
electrons  are  trapped  in  the  ponderomotive  potential  well  or 
"bucket"  produced  by  the  combined  electromagnetic  fields  of 
the  wiggler  and  pulse.  At  this  stage,  the  pulse  continues  to 
grow  with  enhanced  efficiency  as  the  bucket  is  decelerated 
down  the  wiggler  and  energy  extracted  from  the  trapped 
electrons.  A  stationary  state  is  reached  when  the  energy 
extracted  from  the  electrons  is  balanced  by  the  energy 
losses . 

Tasks  1,  2  and  3  address  those  issues  of  the  operation 
of  an  FEL  oscillator  which  relate  to  the  startup  phase, 
saturation  of  linear  gain  with  effective  electron  trapping, 
and  stable  propagation  of  the  finite  amplitude  EM  pulse  in 
steady  state.  The  principal  results  are: 

1)  Adequate  linear  gain  can  be  obtained  in  a  variable 
parameter  FEL  oscillator  sufficient  to  overcome  the  losses  at 
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the  mirrors  and  to  grow  the  EM  pulse  to  a  large  amplitude  in 
a  finite  number  of  passes. 

2)  The  linear  phase  of  growth  saturates  at  a  high 
enough  level  to  ensure  trapping  except  for  systems  which  seek 
to  enhance  linear  gain  by  use  of  a  long  constant  parameter 
section  at  the  front. 

3)  At  high  saturation  amplitudes,  frequency 
discrimination  is  required  to  prevent  nonlinear  breakup  of 
the  EM  pulse  due  to  the  coupling  of  unstable  sideband  modes 
to  the  periodic  motion  of  the  trapped  electrons  in  the 
bucket.  Qualitative  agreement  with  the  theory  of  sideband 
instabilities  is  obtained  for  the  width  of  the  frequency 
filter  required  to  suppress  sideband  growth. 

4)  An  FEL  oscillator  with  an  EM  pulse  of  one  micron 
wavelength  has  been  simulated.  With  proper  frequency 
discrimination,  the  EM  pulse  was  grown  from  noise  to  a  stable 
steady  state  and  30  percent  efficiency  was  achieved. 

5)  A  criterion  has  been  derived  (and  verified 
numerically)  for  the  minimum  current  required  to  have 
adequate  linear  gain  as  well  as  frequency  discrimination: 

Ay  r  2  y  2 

<I>  >  1.42  x  105  (1  -  r)  ‘T|2a  2f  r  amps, 

w 

where  the  electron  mici.y  Ase  was  assumed  to  be  ultra-short, 
is  the  current  averted  ove  a  slippage  distance,  rQ 
the  beam  radius,  Ayr  the  change  in  the  resonant  energy  Yr 
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while  traversing  the  wiggler  of  length  L,  aw  the 
dimensionless  vector  potential  amplitude  of  the  wiggler, 
f  the  filling  factor,  and  (1-r)  the  fractional  reduction  in 
amplitude  due  to  energy  losses  on  reflection  at  the  mirrors. 

6)  The  limitation  on  beam  thermal  spread  is  determined 
more  by  the  need  for  good  linear  gain  than  by  that  for 
effective  electron  trapping. 

These  results  were  first  obtained  in  the  limit  of  an 
ultra-short  electron  micropulse  modeled  as  a  6-function. 

Due  to  computer  time  limitations,  this  limit  is  easier  to 
simulate  since  the  number  of  particle  orbits  which  must 
be  followed  is  far  fewer  than  that  required  to  represent  a 
long  continuous  electron  micropulse.  This  limit  is  also 
analytically  tractable.  However,  the  basic  physics  of  the 
electron-photon  interaction  is  qualitatively  the  same  for 
long  electron  micropulses.  The  details  of  this  investigation 
may  be  found  in  the  first  annual  report. 

The  theory  of  linear  gain  for  the  more  realistic  limit 
of  long  electron  micropulses  is  discussed  in  Appendix  A. 
Simulations  of  micropulses  several  slippage  distances  long 
exhibit  a  behavior  similar  to  that  of  ultra-short  pulses. 
Because  of  computer  time  constraints,  only  a  limited  number 
of  simulations  of  finite  length  micropulses  have  been  done 
(up  to  four  slippage  distances  long).  It  was  again  found 
that  frequency  discrimination  is  necessary  to  suppress  the 
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growth  of  unstable  sidebands.  In  the  absence  of  frequency 

discrimination,  the  EM  pulse  shape  at  saturation  is  broad  but 

7  .... 

very  irregular.  With  frequency  discrimination,  smooth  pulse 
shapes  were  obtained,  although  the  saturated  pulse  energy  was 
somewhat  smaller.  However,  the  indications  are  that  with 
longer  micropulses  such  as  those  actually  envisaged, 
frequency  discrimination  would  be  effective  not  only  in 
producing  smoother  saturated  pulse  shapes,  but  also 
relatively  higher  pulse  energies. 

The  physics  of  the  saturated  state  should  not  depend  on 
the  electron  micropulse  length  since  the  EM  pulse  is  many 
synchrotron  bounce  periods  in  length.  From  arguments  similar 
to  those  used  in  the  case  of  the  ultra-short  pulse,  a 
corresponding  minimum  current  can  be  derived  whic  insures 
that  a  value  of  frequency  discrimination  may  be  chosen  to 
allow  for  both  linear  gain  and  stable  propagation: 


2.77  x  104  (1  -  r)  (Ayr) 


Ip  is  the  peak  current.  Note  that  the  usual  resonance 
condition  kg  =  2kwYr2/(l  +  aj  )  applies  and  that  usually 
the  length  L  is  limited  by  the  Rayleigh  condition 


L  <  k  r  2  . 
sro 

Appendix  A. 


Details  of  the  derivation  are  given  in 
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FEL  AMPLIFIER 


In  high-gain  FEL  amplifiers,  electrons  are  trapped  in 
the  bucket  at  the  front  of  the  wiggler,  the  bucket  is 
decelerated  and  energy  is  transferred  from  the  trapped 
electrons  to  the  EM  pulse  in  one  pass.  Since  the  trapped 
electrons  undergo  many  bounce  oscillations  in  the  bucket 
during  its  passage  through  the  wiggler,  the  growth  of  the 
unstable  sideband  modes  which  can  detrap  electrons  and 
reduce  energy  extraction  efficiency  is  also  a  major  concern. 

In  Task  4,  a  series  of  simulation  runs  were  carried 
out  using  the  one-dimensional  long  pulse  simulation  code 
described  in  Appendix  B  in  order  to  assess  the  effect  of 
sideband  instabilities  on  efficiency  in  the  ATA  FEL  amplifier 
experiment  planned  at  the  Lawrence  Livermore  National 
Laboratory  (LLNL).  A  typical  set  of  design  parameters  for 
the  proposed  ATA  FEL  amplifier  is: 

Table  1 

Y  =  98.85  +  1.87 

Ip  =  5.93  kA 

r  =  0.387  cm 

o 

Pinput  =  2,4  Gigawatts  Xw  =  8  cm 

tp  =  0.447  cm  AQw  =  3.026  kilogauss-cm 

As  =  10.6  x  10"4  cm  Ayr/yr  =  0.33 

L  =  22  meters 


down  the  wiggler.  The 


Ayr  is  the  effective  change  in  yr 
actual  wiggler  profile  is  shown  in  Figure  1. 

The  results  of  the  simulations  using  the  above 
parameters  are: 

1.  With  zero  initial  noise  perturbation  in  the 
EM  pulse,  an  FEL  efficiency  of  18  percent  was  observed, 
agreeing  with  LLNL  estimates.  About  54  percent  of  the 
electrons  were  effectively  trapped.  Figure  2  displays 
a  plot  of  the  dimensionless  pulse  amplitude  as  a 
function  of  the  dimensionless  independent  variable 

v  =  (t  -  z/c)  c/L( c/v-1 )  for  u  =  (z/v  -  t)  c/L(c/v-l)  =  0. 
Note  that  kw(z)  =  constant.  The  wiggler  extends  from  v  =  0 
to  v  =  1 .  The  pulse  amplitude  increases  by  a  factor  of  4.8 
through  the  wiggler.  Outside  the  wiggler  (v  >  1),  the  pulse 
amplitude  is  constant. 

The  bounce  frequency  wB  of  electron  trapped  at  the 
bottom  of  the  bucket  varied  from  ^Bi^w^/^s  =  at  the 

front  of  the  wiggler  to  wBf^wL^Ws  =  at 

the  wiggler. 

2.  With  initial  perturbations  of  sideband  modes  at  a 
frequency  w  in  the  range  of  the  bounce  frequency  w 

O 

large  sideband  amplifications  were  observed.  Figure  3 
displays  a  plot  of  the  variation  of  pulse  amplitude  with  v 
for  the  case  where  the  initial  amplitude  of  the  sideband  mode 
at  frequency  lok^L/uig  =  90  was  0.001  of  the  main  pulse 


7.0 


^ - 1 - 1 - 1 - 

-0.3  0.5  1.3  2.1  2.9 

£  t  -  z/V 
V  L  c/V  -  1 

Figure  3.  Normalized  output  optical  pulse  amplitude 

profile  for  a  simulated  FEL  amplifier  - 

10"3  relative  initial  sideband  perturbation 

at  ;juk  L/jj  =90. 
w  s 
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amplitude.  The  sideband  mode  grew  in  amplitude  by  a  factor 
of  ~2,660.  In  Table  2,  the  amplitude  amplification  factor  of 
sideband  modes  is  tabulated  for  different  values  of  frequency 

wk  L/oj  . 

w  s 

Table  2 


Sideband 
Frequency  oj 

k  L  u)  /k  c 
w  s 

Ratio  of  Initial 
Sideband  Amplitude 
to  Main  Pulse  Amplitude 

Ampli f ication 
Factor 

50 

<n 

l 

o 

r-i 

3.7 

X 

102 

90 

10-3 

26.6 

X 

102 

130 

10-3 

4.6 

X 

102 

150 

10‘3 

7.7 

X 

102 

165 

10"3 

5.9 

X 

102 

180 

10-3 

4.5 

X 

102 

The  observed  sideband  amplification  factors  are 
consistent  with  a  theoretical  analysis  of  sideband  growth  for 
a  high  extraction  FEL  amplifier  in  which  the  electrons  are 
all  deeply  trapped.  The  details  of  this  analysis  will  be 
reported  elsewhere. 

3.  With  an  initial  white  noise  perturbation  containing 
0.5%  of  the  optical  pul3e  energy  but  only  about  10-4  of  the 
optical  pulse  energy  in  the  sideband  frequency  range,  the 
sidebands  grew  to  appreciable  amplitudes  and  the  resulting 
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jutput  pulse  shape  (v  >  1)  displayed  in  Figure  4  was  highly 
irregular.  The  efficiency  was  reduced  from  18%  to  15%  and 
the  effective  trapping  from  54%  to  44%. 

No  serious  degradation  in  efficiency  occurs  until  the 
noise  on  the  incoming  laser  or  e-beam  pulse  is  of  order  10-4, 
well  above  that  due  to  spontaneous  emission.  Clearly  a  2-D 
sideband  theory  is  required  to  assess  whether  beam  quality 
under  circumstances  such  as  shown  in  Figure  4  is  severely 
degraded . 

Since  the  growth  of  sidebands  in  an  amplifier  is 
exponential,  the  growth  may  be  expected  to  be  a  sensitive 
function  of  detailed  design  parameters. 

While  even  the  large  levels  of  growth  predicted  here 
would  not  be  sufficient  to  amplify  spontaneous  emission,  it 
is  clear  that  it  is  essential  that  noise  levels  on  the 
electron  beam  or  laser  not  be  excessive.  We  may  note  that 
coherent  transverse  oscillations  of  the  beam  induced  by 
accelerator  or  transport  systems  could  represent  such  a  noise 
source.  The  second  and  higher  harmonics  of  the  betatron 
frequency  are  in  the  dangerous  range  for  sideband  growth. 
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Figure  4.  Normalized  output  optical  pulse  amplitude  profile 
for  a  simulated  FEL  amplifier  -  10-4  relative 
initial  sideband  energy  with  uniform  noise 
spectrum. 
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SIDEBANDS 


One  of  the  principal  problems  with  tapered  wiggler 
oscillators  is  their  instability  with  respect  to  sideband 
growth  necessitating  frequency  discriminatory  optics.  The 
previous  theory  has  neglected  the  effect  of  the  electrons 
on  the  optical  phase.  The  theory  has  now  been  redone 
exactly.  Some  details  are  presented  in  Appendix  C. 

In  Figures  5  and  6,  the  results  of  analytic  theory 
and  simulations  for  linear  growth  rates  we  have  done  are 
displayed.  Growth  of  sideband  minus  growth  of  signal  per 
pass  is  plotted  in  dimensionless  units,  such  that  signal 
growth  in  these  units  should  be  sin  4/  (where  sin  vr 
measures  the  acceleration  of  the  bucket  and  is  defined  below 
in  terms  of  the  Hamiltonian),  versus  sideband  frequency  <  in 
units  of  (a 0>  z ,  the  synchrotron  frequency  of  an  electron  at 
the  bottom  of  the  ponderomotive  well.  The  growth  rates  are 
presented  for  two  particle  distribution  functions. 

(1)  Uniformly  filled  bucket.  In  Figure  5,  linear  gain 
curves  for  the  old  and  new  theories  are  compared. 

While  Case  (1)  is  easy  to  calculate  analytically,  it  is 
impossible  to  simulate  because  of  nonlinear  effects  on 
particles  at  the  separatrix.  Thus  we  look  also  at  a 
distribution  f  with  a  linear  gradient. 
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(2)  f0  =  (Hmax  -  H)(Hmax  ~  Hmin>'  where  H  is  the 
particle  Hamiltonian  H  =  -a  (cosiJj  +  tj;sin  'l>r )  .  The  solid 

curve  in  Figure  6  is  the  numerically  evaluated  analytic 
dispersion  relation  for  this  case  with  sin  <Pr  =  0.  it  will 
be  observed  that  for  sin  tj^  =  0,  the  curves  for  case  (1)  and 
case  (2)  are  quite  similar  so  that  results  may  be  expected  to 
be  prototypical  for  all  distributions  which  extend  over  the 
whole  bucket.  Sideband  growths  would  be  larger  for  distribu¬ 
tions  concentrated  near  the  bottom  of  the  bucket. 

The  points  are  simulation  results  for  sin  =  o,  0.2, 
0.5.  It  will  be  observed  that  for  sin  y  =0.2  and  0.5,  the 
peak  growth  of  sidebands  is  4.5  and  2.5,  respectively,  times 
signal  growth.  More  detailed  comparison  between  theory  and 
simulation  for  sin  ip^  ~  o  is  given  in  Appendix  C. 

It  should  be  noted  that  for  sin  -pr  ?  o,  the  new 
theory  predicts  slightly  higher  growth  than  that  given 
in  Reference  2.  In  particular,  sideband  growth  exceeds 

signal  growth  at  frequencies  close  to  that  of  the  signal, 

2 

5=1+  &(</n)  .  Thus,  some  frequency  discrimination  is 
required  even  against  low-lying  sidebands. 
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IV.  GAIN-EXPANDED  FELs 


In  a  gain-expanded  FEL,  the  sensitivity  of  gain  to 
variations  in  the  energy  of  the  exciting  electrons  is  reduced 
by  using  a  periodic  magnetic  wiggler  with  transverse 
gradients  such  that  the  electrons  which  lose  energy  move  to 
different  transverse  positions  and  continue  to  remain  in 
resonance  with  the  ponderomoti ve  potential  well.  The  reduced 
sensitivity  to  energy  makes  it  attractive  to  operate  a  gain- 
expanded  FEL  in  conjunction  with  a  storage  ring  since  many 
passes  of  the  circulating  electrons  through  the  FEL  would 
occur  before  gain  is  significantly  affected  by  energy  spread. 
However,  transverse  betatron  oscillations  of  the  electron 
trajectories  are  simultaneously  excited  and  may  negate  the 
effectiveness  of  gain  expansion  unless  the  oscillation 
amplitudes  are  kept  small.  For  steady  state  operation  in  a 
storage  ring,  it  is  desirable  to  maximize  the  gain  while 
minimizing  the  level  of  betatron  oscillations. 

In  order  to  obtain  some  insight  into  the  gain 
characteristics  of  gain-expanded  FELs,  a  Hamiltonian  formula¬ 
tion  of  the  equations  of  motion  for  the  "thin  lens" 
gain-expanded  FEL  was  derived.  Two  fundamental  limitations 
on  gain-expansion  schemes  were  established:  The  Manley-Rowe 
relation  which  relates  gain  to  betatron  excitation,  and  the 
generalized  Madey  theorem  relating  gain  to  energy  spread  and 
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betatron  excitation.  The  details  of  this  analysis  can  be 
found  in  Appendix  A  of  the  1982  Annual  Report. 


These  analytic  results  were  compared  this  year  with 

Q 

John  Madey's  Monte  Carlo  simulation  of  gain-expanded  FELs  in 
storage  rings  and  the  following  conclusions  were  reached: 

1.  The  latest  version  of  the  simulations  are  now 
consistent  with  the  Manley-Rowe  and  the  gain-spread 
relations . 

2.  At  low  laser  power,  there  is  good  agreement  between 
the  simulation  and  the  analytic  results.  At  higher  power, 
the  Monte  Carlo  simulations  exhibit  some  degradation  in 
performance.  This  may  be  due  to  cavity  coupling  to  betatron 
oscillations  which  lead  to  resonance  spreading,  an  effect  not 
included  in  the  analysis.  However,  proper  magnet  design 
which  minimizes  this  coupling  could  improve  the  performance. 

3.  While  all  designs  are  limited  in  consequence  of  the 

above  theorems  to  AE,  /Ae  ,  <  oy/T  where  77 

laser'  ^synchrotron  ]/  ' 

is  the  rms  spread  in  electron  energy,  it  would  appear  that 
gain-expanded  wigglers  may  be  advantageous  compared  to 
conventional  wigglers  by  virtue  of  larger  energy  acceptance. 
Studies  of  various  designs  will  continue  in  collaboration 
with  Madey's  group. 
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V.  PHASE  AREA  DISPLACEMENT 


In  the  deceleration  of  relativistic  beam  electrons  by  a 
phase  area  displacement  wiggler,  the  resonant  energy  of  the 
wiggler  is  increased  from  the  front  to  the  back  such  that  the 
bucket  is  moved  through  the  phase  area  occupied  by  the 
electrons.  The  result  is  a  downward  displacement  in  the 
electron  energy  by  an  amount  of  the  order  of  the  phase  area 
of  the  bucket  divided  by  2t.  In  the  idealized  limit  of 
infinitely  slow  deceleration,  there  should  be  no  energy 
spreading . 

This  method  of  beam  energy  extraction  has  two  features 
which  are  attractive  with  respect  to  operating  an  PEL 
oscillator  in  conjunction  with  a  storage  ring:  The  energy 
extracted  is  insensitive  to  the  beam  energy  spread,  and  the 
ratio  of  the  energy  extracted  to  the  increase  in  root  mean 
square  energy  spread  can  be  made  to  be  small. 

Numerical  simulations  of  a  one-dimensional  FEL  with  a 
phase  area  displacement  wiggler  have  been  done  and  the 
results  were  consistent  with  the  theoretical  picture  of  an 
average  energy  loss  and  an  average  energy  spread  independent 
of  the  initial  energy,  although  the  simulated  energy  loss  was 
about  20  percent  lower  and  the  spread  about  25  percent  higher 
than  predicted  by  the  crude  theory.  The  details  of  these 
simulations  are  described  in  Appendix  D. 
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An  FEL  oscillator  can  be  operated  in  steady  state  with 
a  storage  ring  if  the  laser-induced  increase  in  energy  spread 


can  be  balanced  by  a  decrease  due  to  incoherent  synchrotron 
radiation  and  the  beam  energy  boosted  by  a  radio- frequency 
cavity  to  compensate  for  the  losses  in  the  wiggler  and  in 
synchrotron  radiation. 

In  order  to  obtain  good  efficiency  during  steady  state 
operation,  it  is  essential  that  the  electron  interaction  with 
the  bucket  be  adiabatic  and  that  electron  trapping  be 
negligible.  Typically,  adiabaticity  requires  long  wigglers. 
The  bucket  must  remain  essentially  constant  through  the 
wiggler  since  bucket  variations  can  affect  not  only 
adiabaticity,  but  lead  to  electron  trapping.  Thua.  it  wi l ' 
be  necessary  to  have  a  smooth  and  long  (many  tines  the 
slippage  distance)  EM  pulse,  entailing  a  good  frequency 
discrimination. 

If  these  conditions  are  satisfied,  the  FEL  efficiency 
(measured  in  terms  of  the  ratio  of  the  energy  extracted  from 
the  electrons  ^Yiaser  to  the  energy  loss  in  synchrotron 
radiation  ^YSyn)  is  estimated  to  be: 


A  Y  _ 

laser 

Ay 

syn 


~  0.05  k  2L2 
w 


V; 


\ 


23 


I 


Here  ~  l  is  the  wiggler  dimensionless  magnetic  vector 

potential,  the  wave  number,  ks  the  optical  wave  number, 

and  r^  the  optical  pulse  radius. 

For  a  particular  design  discussed  in  Appendix  D,  a 

value  of  ^^laser^ ^Ysyn  ~  0.15  was  observed  in  the 
simulation,  in  rough  agreement  with  theoretical  predictions. 

Long  wigglers  containing  many  wiggler  periods  and  EM 
pulses  with  large  circulating  peak  powers  are  particularly 
effective  in  obtaining  high  efficiencies. 

The  linear  gain  per  pass  tends  to  be  very  small  unless 
appreciable  micropulse  peak  current  densities  are  available. 
Detailed  formulae  are  given  in  Appendix  D.  Even  if  linear 
gain  is  adequate,  it  is  not  clear  whether  the  desired  final 
steady  state  is  accessible  by  growing  the  EM  pulse  from  noise 
levels.  The  scenario  of  the  time  evolution  to  a  steady  state 
is  complex  and  remains  to  be  elucidated.  It  is  thus  not  yet 
clear  whether  a  phase  displacement  pulse  can  be  grown  from 
noise  or  whether  some  alternate  startup  strategy  must  be 


found . 


VI.  TWO-DIMENSIONAL  EFFECTS 


Hitherto,  the  electron  equations  of  motion  have  been 
considered  as  a  one-dimensional  problem  in  which  the  electron 
motion  is  reduced  to  the  pendulum  equation  in  the  ponderomo- 
tive  potential  well  formed  by  the  wiggler  and  EM  pulse 
fields . 

In  fact,  the  electromagnetic  fields  do  depend  on 
transverse  dimensions.  Transverse  variations  of  the  wiggler 
field  produce  transverse  betatron  oscillations  and  diffrac¬ 
tion  effects  introduce  curvature  of  the  wave  front. 

In  the  simplest  approximation,  the  transverse  betatron 
oscillation  ("frequency"  kg)  is  decoupled  from  the 
longitudinal  "synchrotron"  oscillation  (bounce  "frequency"  H) 
in  the  ponderomoti ve  potential  well.  However,  curvature  of 
the  wave  front  couples  the  two  oscillations.  This  coupling 
is  weak  except  near  resonance  when  2k^  =  i!  which  occurs 
typically  at  peak  circulating  power  levels  near  1  Gigawatt. 
The  details  of  the  effect  of  this  resonance  coupling  in 
detrapping  electrons  is  described  in  Appendix  E. 

Under  resonant  conditions,  the  coupling  is  marginally 
strong  enough  to  lead  to  some  detrapping.  Thus,  a  detailed 
numerical  simulation  of  specific  cases  may  be  required. 


In  an  amplifier,  it  appears  that  the  parameters  may 
change  rapidly  enough  that  the  resonance  is  passed  through 


without  significant  detrapping.  During  the  buildup  phase  of 
an  oscillator  passing  through  resonant  power  levels,  the  peak 
potential  detrapping  is  estimated  to  be  of  the  order  of 
40  percent  if  the  emittance  is  the  maximum  allowed  by  other 
considerations.  Hence,  if  a  reasonable  gain  margin  exists, 
this  detrapping  is  not  essential.  With  smaller  emittances, 
the  effect  is,  of  course,  smaller.  For  actual  cylindrical 
cases  where  many  orbits  do  not  pass  close  to  the  axis,  the 
effect  is  further  reduced. 

In  the  analysis,  the  self-consistent  effect  of  betatron 
motion  on  emission  was  not  studied,  but  only  the  effect  of  an 
assumed  wave  shape  on  electron  trapping.  Subject  to  these 
caveats,  it  would  appears  that  two-dimensional  effects  do  not 
seriously  perturb  the  simple  one-dimensional  picture  of  an 
FEL  or  introduce  significant  additional  design  constraints. 


ACKNOWLEDGEMENTS 


We  would  like  to  acknowledge  helpful  discussions  with 
Norman  Kroll  and  Lee  Sloan  concerning  phase  area  displacement 
FEL  operation,  and  with  Norman  Kroll  concerning  sideband 
instability  studies.  We  would  also  like  to  acknowledge 
useful  discussions  with  C.  Brau,  W.  B.  Colson,  J.  C. 
Goldstein,  G.  Neil,  j.  Madey,  D.  Prosnitz,  J.  Slater, 

P.  Sprangle  and  C.  M.  Tang. 


27 


I 


REFERENCES 


1.  L.  R.  Elias,  W.  M.  Fairbank,  J.  M.  J.  Madey,  H.  A.  Schwettman 
and  T.  I.  Smith,  Phys.  Rev.  Lett.  35_,  pp.  717-720  (1976). 

B.  E.  Newman,  R.  W.  Warren,  K.  Boyer,  J.  C.  Goldstein,  M.  C. 
Whitehead  and  C.  A.  Brau,  "Optical  Gain  Results  of  the  Los 
Alamos  Free  Electron  Laser  Amplifier  Experiment,"  Submitted 
to  Society  of  Photo-Optical  Instrumentation  Engineers, 

Special  Publication  No.  453  (1983  Free  Electron  Laser 
Workshop  -  Orcas  Island,  WA) . 

2.  N.  M.  Kroll,  P.  Morton  and  M.  N.  Rosenbluth,  "Free  Electron 
Lasers  with  Variable  Parameter  Wigglers,"  IEEE  Journal  of 
Quantum  Electronics,  Vol.  QE-17,  pp.  1436-1468  (1981). 

3.  T.  I.  Smith,  J.  M.  J.  Madey,  L.  R.  Elias  and  D.  A.  G. 

Deacon,  J.  Appl .  Phys.  j^O,  4580  (1979). 

N.  M.  Kroll,  P.  L.  Morton,  M.  N.  Rosenbluth,  J.  N.  Eckstein 
and  J.  M.  J.  Madey,  "Theory  of  the  Transverse  Gradient 
Wiggler,"  IEEE  Journal  of  Quantum  Electronics,  Vol.  QE-17, 
pp.  1496-1507  (1981). 

4.  M.  N.  Rosenbluth,  H.  Vernon  Wong  and  B.  N.  Moore,  "Annual 
Technical  Report  for  Theoretical  Studies  on  Free  Electron 
Lasers,"  ARA  Report  No.  I-ARA-82-U-89 ,  August  1982. 

M.  N.  Rosenbluth,  H.  Vernon  Wong  and  B.  N.  Moore,  "Free 
Electron  Laser  (Oscillator)  -  Linear  Gain  and  Stable 
Pulse  Propagation,"  Submitted  to  Society  of  Photo-Optical 
Instrumentation  Engineers,  Special  Publication  No.  453 
(1983  Free  Electron  Laser  Workshop  -  Orcas  Island,  WA) . 

5.  D.  Prosnitz,  A.  Szoke  and  V.  K.  Neil,  "One-Dimensional 
Computer  Simulation  of  the  Variable  Wiggler  Free  Electron 
Laser,"  Vol.  7,  Free  Electron  Generators  of  Coherent 
Radiation,  Addison-Wesley,  Reading,  MA,  pp.  571-587  (1980). 

D.  Prosnitz,  private  communication. 

6.  M.  N.  Rosenbluth  and  H.  Vernon  Wong,  "Annual  Technical 
Report  for  Theoretical  Studies  on  Free  Electron  Lasers,” 
Appendix  A,  ARA  Report  No.  I-ARA-82-U-89,  August  1982. 

7.  John  C.  Goldstein  and  W.  B.  Colson,  "Control  of  Optical 
Pulse  Modulation  Due  to  the  Sideband  Instability  in  Free 
Electron  Lasers,"  Submitted  to  the  Proceedings  of  the 
International  Conference  on  Lasers  1982,  New  Orleans,  LA, 
December  15-17,  1982. 

8.  J.  Madey,  private  communication. 


28 


ft 


APPENDIX  A 


THEORY  OF  LINEAR  GAIN:  FREE  ELECTRON  LASERS  OPERATED 
IN  OSCILLATOR  MODE  -  FINITE  LENGTH  ELECTRON  PULSE 

In  the  1982  Annual  Report,1  the  theory  of  linear  gain 

was  investigated  for  a  Free  Electron  Laser  (FEL)  operated 

as  an  oscillator  with  ultra-short  electron  pulses.  In  such 

a  device,  a  periodic  magnetic  field  wiggler  of  length  L 

and  wave  number  k  is  positioned  between  mirrors  so  that 

w 

an  electromagnetic  pulse  can  be  reflected  repeatedly 

through  the  wiggler.  The  FEL  is  driven  by  a  successive 

series  of  ultra-short  electron  beamlets  injected  at 

periodic  intervals  so  that  there  is  overlap  of  beamlet  and 

the  electromagnetic  pulse  on  each  forward  pass  through 

the  wiggler.  By  ultra-short  is  meant  beamlets  with 

length  l  much  less  than  the  slippage  distance  between 

electrons  and  EM  pulse  s  *  kwL/ks ,  with  kg  the  wave 

number  of  the  EM  pulse.  In  fact  as  we  will  see  herein  for 

the  results  of  Reference  1  to  be  valid,  we  must  require 
u 

s/T  >  l  >  2-rr/k  with  T  =  2k  L  Ay  /y  >>  1,  where 

'  s  w  r  r 

£y r/ is  the  fractional  change  in  the  resonant  energy  down 
the  wiggler.  The  linear  eigenmode  equations,  obtained  by 
linearizing  the  FEL.  equations  in  the  amplitude  of  the  laser 
pulse,  were  analyzed  and  the  linear  gain  per  pass  determined 
as  a  function  of  the  FEL  parameters. 
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In  this  paper,  we  extend  the  analysis  to  cover  the 

case  of  electron  beamlets  with  lengths  much  larger  than  the 

t 

slippage  distance,  l  >  s/T*.  our  starting  point  will  be 
the  conventional  set  of  FEL  equations  applicable  in  the 
Compton  regime,  as  derived  and  discussed  in  Reference  1, 
These  will  be  restated  without  derivation  in  Section  1  for 
completeness  of  presentation. 

We  anticipate  the  need  to  suppress  the  growth  of 
sideband  instabilities  when  the  pulse  reaches  a  finite 
amplitude  sufficient  to  trap  beam  electrons.  This  may  be 
accomplished  by  frequency  discrimination  of  the  pulse  in 
which  frequencies  above  and  below  the  desired  pulse 
frequency  ojs  are  attenuated.  Thus,  the  analysis  will 
encompass  variable  parameter  wigglers  without  and  with 
frequency  discrimination.  The  linear  eigenmode  equation 
and  their  solutions  are  discussed  in  Section  2  and 
Section  3  for  wigglers  without  and  with  frequency  discrimi¬ 
nation,  respectively. 

In  Section  4,  the  results  are  compared  with  those 
previously  obtained  for  the  ultra-short  electron  beamlets. 

I .  FEL  Equations 

The  set  of  equations  which  describe  the  temporal 
(t)  and  spatial  (z)  evolution  of  the  electromagnetic  pulse 
in  the  wiggler  are  most  conveniently  expressed  in  terms  of 
the  new  independent  variables  u  and  v: 


where  V  =  w  /  (5<w+^s )  is  the  electron  resonant  longitudinal 
velocity . 

The  periodic  magnetic  field  wiggler  and  the 
electromagnetic  pulse  are  represented  by  circularly 
polarized  vector  potentials  A^  and  Ag  with  negligi¬ 
ble  spatial  variations  transverse  to  the  direction  of 
electron  beam  propagation.  The  case  of  plane  polarization 
is  given  by  letting  A  -*•  h/JY.  The  electromagnetic  pulse  is 
approximated  by  a  plane  wave  ~Ag  exp(-iws  t  +  ikg  z  +  i£) 
with  slowly  varying  amplitude  Ag  and  phase 
The  electron  equations  of  motion  are: 

=  v  (1) 

3v 

ii  ,  r  ♦  Uiae1*  -  (2) 

3v  2  [_  J 
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wh  ere 


2k  L 

-  w 

_  -  ( >  —  Y  ) 

Y_  r 


ip  =  fkdz+kz-wt 
J  w  s  s 

0 

v&f  (l 


w 


2k  k  L2 

-  w  s  iC 

a  *  - ; -  a  a  e 

y  2  w  s 


2k  L 
r  =  -J±- 


y  and  ip  are  the  electron  energy  and  phase  angle. 

a  =  eA  /me2  and  a  =  eA  /me 2  are  the  dimensionless 
w  w  s  s 

vector  potential  amplitudes  of  the  magnetic  field  wiggler 

and  pulse,  respectively.  a*  is  the  complex  conjugate  of 

a.  In  these  equations,  a  and  k  are  considered  to 

w  w 

be  constants,  except  as  they  determine  the  resonant  energy 

Yr ,  and  the  variable  parameter  wiggler  is  modeled  by 

r  jt  0,  where  uy^  corresponds  to  the  change  in  the  resonant 

energy  in  going  through  the  wiggler.  We  are  engaged 

primarily  with  high  efficiency  cases,  T  >>  1 .  It  is  also 

assumed  that  y  -  y  < <  y 

r  r  • 

The  equations  for  a^  and  ;  (from  Maxwell's 
equations)  are: 
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:  a 

3  U 


(3) 


where 


8Tre2  N  k  L2  a2 
_ T  w  w 


Y  me  V  ( 1  +  a  2  ) 
r  w 


/ 


du  h(u) 


1 


Nt  is  the  total  number  of  electrons  in  the  beamlet  per 
unit  area,  N^,  =  J dz  n(z,t)  where  n(z,t)  is  the  beam 
density.  h(u)  is  a  form  factor  which  is  determined  by  the 
profile  of  the  electron  beamlet.  For  an  ultra-short  pulse 
as  described  in  Reference  1,  h(u)  =  6(u).  The  angular 
brackets  imply  integration  over  the  initial  energy  distri¬ 
bution  and  average  over  the  random  optical  phase  of  the 
electrons 


V 


where  f(y  )  is  the 

o 

electrons . 

In  u-v  space, 
the  lines  u  +  v  =  0 


initial  energy  distribution  of  the 

the  wiggler  (L  2  z  20)  lies  between 
and  u  +  v  =  1.  The  beam  electrons 
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move  on  lines  of  constant  u  and  the  photons  of  the 
electromagnetic  pulse,  propagating  in  the  beam  direction, 
move  on  lines  of  constant  v.  The  beam  electron  and 


photons  interact  only  when  their  trajectories  in  the  u,v 
plane  intersect  within  the  lines  u  +  v  =  0  and  u  +  v  =  1. 

II .  Linear  Theory  -  Without  Frequency  Discrimination 

The  formulation  of  the  linear  eigenmode  equations 
has  previously  been  discussed  in  Reference  1.  If  we  follow 
the  same  procedure,  we  obtain  for  the  pulse  amplitude  a11 
on  the  nth  pass: 

_  „  n 

~  n  _  3a 

(1  -  r  +  5)  a  -  6  — 


1-v 


=  in 


duh(u)  <^exp(-iCO^> 


(4 


-v 


<^exp  (-iv))> 


1 

2 


f  ( y  ) 
o 


(v-v ’ ) 


-u 


-  n 
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(V  ) 


(  5 


exp 


(v’-v)  +  i 


(v'+u)‘ 


-  i 


_  (u+v)‘ 


34 


l 


ft 


r  accounts  for  the  reduction  in  amplitude  due  to  er.erov 
losses  on  reflection  at  the  mirrors  and  is  close  to  but 


less  than  unity.  £  =  k  CAt/k  L  <<  1  where  At  is  the 

s  w 

pass-to-pass  temporal  advance  of  the  electromagnetic  pulse 
relative  to  the  electron  beamlet  on  entry  into  the  wiggler, 
and  yq  is  the  value  of  y  at  wiggler  input.  We  have 
expanded  a(v+3)  =  a  +  8da/dv  +  82  .  .  .  .  The  linear  gain 

and  phase  shift  per  pass  6  is  the  eigenvalue.  The 
boundary  conditions  are  an(v)  -*  0,  v  ■+  t  °°. 

In  deriving  Equation  (4)  and  Equation  (5),  it  is 
assumed  that  <5  <<  1. 

Substituting  Equation  (5)  in  Equation  (4) 


od 

OV 


-  (1  -  r  +  5)  a 


h  (y-v)  a  (v-x) 


I  (x)  x 


where 


exp 


iTxv 


nh (-v) 
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f dx  a (v-x)  I (x) 


i .  x 


exp 
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.  X  ; 


)  exp 
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and  h(y-v)  is  assumed  in  accordance  with 
assumption  to  be  a  slowly  varying  function 
so  that  h(v-v)  ss  h(-v)  +  1  >  y  >  0. 

the  superscript  n  on  a  . 

If  we  take  h(u)  to  be  Lorentzian 

1 


“oj 

and  we  define  the  Fourier  transform  b(k)  of  a: 


1 

h(u)  =  - 

_u 

o 


the  long-pulse 
of  its  argument 
Hereafter  we  drop 


b  ( k) 


+  ■= 


-i  k  v 


qv  e 


a  ( v) 


—  oc 


we  obtain  from  Eauation  (f): 


wh  ere 


A(0 


G  (<) 
1 


6  =  (1  -  r  +  6)/Br!'J 

i 

n 

ai  "  2ttu  er2 

0 


The  parameter  is  a  measure  of  the  length 

of  the  pulse  (full  width  at  half  maximum),  u  =  £/L( 1-V/c) , 

o 

where  L(l-V/c)  is  the  "slippage  distance." 

We  see  from  the  dimensionless  form  of  Equation  (7) 

that  the  long  or  short  pulse  approximations  depend  on 

whether  Fi;u^  is  large  or  small  compared  to  unity.  It  may 

V 

)..<_•  seen  a  posteriori  that  F  ‘u  >1  justifies  our  expansion 

o 

of  h( y-v) . 

If  G  (<)  is  now  expanded  about  <  =  < c  where 

cG  ,  (<  ) 

o 

3  K 


=  0 


(9) 
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Equation  (7)  may  be  approximated  by: 


i 


1  +  G  (<  )  + 

|  ,  0 


(10) 


Furthermore,  if 


1  +  G  (<  ) 
1  c 


0 


(ID 


Equation  (10)  will  have  the  solution: 


and  this  guarantees  a(v)  -*  0  as  v  -*■  t  <»  provided  that 


Re 


> 


0 


(12) 


(  12)  , 


Equation  (9)  and  Equation  (11),  subject  to  Equation 
constitute  the  eigenvalue  equation  which  determines  6  . 
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If  we  take  f(v„)  to  be  Lorentzian 


f  h 


we  obtain  for  K^/f  )'• 


(A)  - 


-e  £ 


5  <Yth Y  iY  >/r" 


and  in  the  limit  of  an  infinitely  long  wiggler  where 
T  >>  1,  Equation  (8)  for  Gi(<)  may  be  approximated  by: 


G  (<) 
i 


(i-1) (i<-:  . ) 


m) 


(13) 


where  Z(6)  is  the  Plasma  Dispersion  Function 


it 


Z  ( r )  =  2i  exp (-6  *  ) 


d  i  exp ( -  - 2 ) 


j 

—  oc 


If  we  substitute  Equation  (13)  for  G;(<)  in 
Equation  (8)  and  Equation  (11),  we  obtain  the  eigenvalue 
eauati ons 


,  1  I 

1  +  rr  z  (r  ) 
2i  o 


a,  Z  <VZ  <V 


8  7u  ‘  Z(-  ) 


=  0 


(14) 
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(i-l)Z(c  ) 

-  U"i)0O  +  z -  C 

0 


(15) 


where 


i  <  +  £ 
0 

*0  "  1  -  i 


Re 


-  a 


»i  I 

Z  Z 


1  z 


>  0 


Eauation  (14)  is  solved  for  6  ,  and  substitution  of  e 

‘  o  o 

in  Equation  (15)  yields  6  as  a  function  of  a^ , 

Fu  2  and  e . 

In  the  limit  of  Fu  2  -  <*>  and  a  E  h/2ttu  £F2  >>  1, 

0  10 

we  use  the  asymptotic  form  of  Z ( S ^ )  in  Equation  (14) 

to  determine  9  2  s  ia  /2,  and  after  substitution  in 

o  i 

Equation  (15),  we  obtain  for  the  linear  gain: 


n-v  +  i)  *  2  (2^7)  ~  ( 


th 


+  1  V 


)  = 


(16) 


More  generally,  Equation  (14)  and  Equation  (15)  are 
solved  numerically  and  the  variation  of  the  dimensionless 
linear  gain 


40 


Re 


5  +e 

i 


=  (1  -  r  +  6) 


2^u  r 


V2 


is  plotted  as  a  function  of  l/ctj  for  various  values  of 

r*u  in  Figure  1.  The  linear  gain  (  y. .  =  0,  r  u  •*  ®) 
o  th  o 

reaches  a  maximum 


Re  l-r+6)  a  1.16  - ! - 

max  _  v, 

TTU 

0 

at  a  value  of 

2ttu  6T2 

- - -  as  3. 

n 

The  curves  for  f^u  =  1,3  were  obtained  by 

o 

numerically  integrating  Equation  (7)  subject  to  the 
boundary  condition  A(<)  -*■  0 ,  <-*■*<». 

There  also  exist  additional  solutions  of 
Equation  (14)  and  Equation  (15)  corresponding  to  eigen- 
modes  with  an  increasing  number  of  nodes  in  the  amplitude 
a(v).  The  linear  gain  for  the  three  "lowest"  eiaenmodes  in 

Vs 

the  limit  of  r  u  >>  1  are  displayed  in  Figure  2.  It  may 

o 

be  seen  that  for  extremely  long  wigglers  these  higher  modes 
may  grow  even  when  the  fundamental  is  damped. 
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With  frequency  discrimination,  the  electromagnetic 
pulse  is  filtered  after  each  pass  through  the  wiggler  by  a 
band  pass  filter  with  frequency  half  width  Ate: 

vc 

Aui  =  - 

L (c/V-1 ) 

where  v  is  an  additional  parameter  which  characterizes 
the  magnitude  of  the  frequency  half  width.  Since  filtering 
leads  to  an  effective  retardation  of  the  pulse  by  l/v 
relative  to  the  electron  beamlet,  while  8  represents  a 
forward  shift,  it  appears  optimal  to  choose  v  =  l/B  to 
avoid  the  pulse  running  ahead  or  lagging  behind  the 
electron  beamlets. 

With  this  choice  of  v  =  1/8,  the  linear  eigenmode 
equation  for  the  nth  pass  is  determined  by: 


-v 


A  straightforward  repetition  of  the  analysis 
described  in  Section  2  yields  the  following  equation  for 
the  Fourier  amplitude  b{<): 
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.  U' 


J  <  * 


B  ( •■• ) 


(  1  *  G-  (*' ))B  (v  ) 


0 


I 

where 


k 


2 


< 


B(<)  =  (<2  +  6  )  £ 

2  2 


r  2 


*«,)  ja{  j(ri) 


exp  (-i  <  £) 


7TU  6 2  r  5/2 

0 


2  (1  -  r  *  6) 
a  z  r 


Thus,  the  eigenvalue  equations  are  given  by 


=  0 


/  1 


G 


2  \ 


1  +  G. (<  )  + 


1 

2 


subject  to  the  requirement 


■•(3) 


>  0 


(22) 


For  an  initial  distribution  function  £(y which 


is  Lorentzian,  and  in  the  limit  of  an  infinitely  long 
wiggler  where  T-  >>  l,  Equation  (19)  for  G2(<)  may  be 


approximated  by: 


G  {<)  =  - 

2 


/ i<  +  £  \ 

M— ) 

( < 2  +  S  2 )  ( i— 1 ) 


(23) 


If  we  substitute  Equation  (23)  for  G„(<)  in 
Equation  (20)  and  Equation  (21),  we  obtain: 


a  2  (5  ) 

2  C 


4  < 


(i-1) a.  Z  (e  ) 

i6ru2  <  z  ( f  ) 
0  0  0 


/ 

0 


=  o 


(24) 


2 (1+i) <„  Z (6  ) 


Z  '  t  ) 
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where 


e 

o 


i<  +  e 

1  -  i 


Re 


( i— l )  z  ’  ( e  ) 

c 


k  Z(3  ) 

o  o 


2k 


o 


+ 


( 1  —  i ) 

4 


>  0 


Equation  (24)  is  solved  for  Bq  ,  and  substitution 

of  6  in  equation  (25)  yields  6  as  a  function  of  a  , 

0  2  2 

.ru2  ,  and  e  . 

c 

In  the  limit  of  -  00 ,  £  -*■  0,  and 

a  =  n  /it  u  62r  *  >>1,  we  use  the  asymptotic  form  of  Z(6  ) 

2  0  0 

in  Equation  (24)  to  determine  S^3  as  a^d-D/8,  and  after 
substitution  in  Equation  (25),  we  obtain  for  the  linear 
gain  <5: 


(1  -  r  +  6) 


(26) 


In  the  limit  of  fu.  —  £  -*  0,  and 

U 

l/Tu  2  <<  1,  we  use  the  series  expansion  of  Z(0  )  in 

0  2  0 

Equation  (24)  and  Equation  (25)  to  determine  6^ ss (1-i) a^/4 
and  the  linear  gain 


(1  -  r  +  5) 


(l-i)n 


(27) 
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Equation  (24)  and  Equation  (25)  are  solved 
numerically,  and  the  variation  of  the  dimensionless  linear 
ga  in 

(  2ttu  r3/2  ) 

Re  ( 6  /a  )  =  Re  (l-r  +  6)  - - - 

2  (  n  ) 

is  plotted  as  a  function  of  l/a22  for  various  values  of 

V '2  W 

7  u  in  Figure  3.  The  linear  gain  (e=0,  Y  u  -+  °°)  reaches 

o  o 

a  maximum 

Re  ( 1  -  r  +  6  )  a:  0.56  - 

mu  r  '2 
o 

for  a  value  of  6  given  by 

1 

2 

a:  C .  81 

k 

The  curves  for  r  u  =  1,3  were  obtained  by 
numerically  integrating  Equation  (18)  subject  to  the 
boundary  condition  B(<)  -»  0 ,  <-*•  +  «». 

In  Figure  4,  the  reduction  in  linear  gain  due  to  a 
finite  energy  spread  y ^  is  displayed,  while  Figure  5 
indicates  that  linear  gain  is  maximized  if  the  input  beam 
energy  is  below  the  resonant  energy. 
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.n  with  Frequency  Discrimination 


Discussion 


by 


The  Iona  electron  beam  pulse  limit  is  characterized 


2k  L 


w  w 


1 


2 


u  >>  1 

o 


where  un  is  the  ratio  of  the  pulse  length  (full  width  at 
half  maximum)  to  the  slippage  distance. 

In  the  limit  of  a  long  wiggler  >>  1)  and  a  long 

electron  beamlet  ( T^u^  >>  1),  the  linear  gain  in  the 
absence  of  frequency  discrimination  is  a  maximum 


Re  ( 1  -  r  +  c  ) 

max 


1.16 


n 


7. 

ttu  r  ‘ 
0 


for  a  value  of  S  given  by  [2^u  SrVo]  ~  3.0,  and  it  is 

c 

2 

close  to  that  estimated  by  Kroll,  et  al .  With  frequency 
discrimination,  the  linear  gain  is  slightly  lower  and  has  a 
max imum 


Re  ( 1  -  r  +  6  ) 

max 


0. 56 


h 

tu  r  2 
0 


for  a  value  of  6  given  by 

.  7-  1  k 


"U  “ 


=  0.85 
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The  corresponding  results  of  the  linear  gain 

k 

analysis  for  the  ultra-short  electron  beamlet  (r  u  <<  1) 


are  Re ( l-r+  o ) 


•  18n/r  at  2  6  r  ^  n  ~  *45  (without 


freauency  discrimination)  and  Re(l-r+6)  ~  .105  n/T 

J  max 

at  6T/n  ‘  ~  .43  (with  frequency  discrimination).  Note  that 

k 

for  pulse  length  F‘  u  ~  1,  the  growth  rates  of  long  and 

c 

short  pulses  are  comparable.  For  longer  pulses  the  growth 

depends  on  peak  current  n/u  • 

o 

Finite  temperature  reduces  gain  significantly  when 

_  ^ 

Y  ,  /T~  >  1.  Linear  gain  may  be  increased  somewhat  (in  the 
th 

presence  of  frequency  discrimination)  if  the  input  beam 
energy  is  decreased  below  the  resonant  energy. 

The  need  for  frequency  discrimination  to  prevent 
nonlinear  breakup  of  the  saturated  large  amplitude  EM 
pulse  due  to  the  growth  of  sideband  instabilities  has 
previously  been  discussed  for  the  case  of  ultra-short 
electron  beamlets.^  It  is  anticipated  that  a  similar 
requirement  will  also  be  necessary  for  the  case  of  long 
electron  beamlets. 

To  suppress  sideband  frequencies  of  the  order  of  a 
quarter  of  the  electron  bounce  frequency  —  ^ | a j  '  ,  the 
reduction  in  the  pulse  energy  due  to  frequency  discrimina¬ 
tion,  that  is, 


i 


should  exceed  the  relative  gain  due  to  sideband 
instabilities.  This  implies 

-S— r  >  2  (1  -  r)  (G  -  1) 

16  v  * 

where  G  is  the  ratio  of  the  sideband  instability  gain  to 
the  EM  pulse  amplitude  gain  at  saturation.  If  we  substitute 
for  the  pulse  amplitude  at  saturation 

i  ~  1 1  n  r 

'a  (l-r)tuQ 

and  we  take  v  =  1/6.  we  obtain 


6 2 


rn 

mu  (1-r)  3 


>5 


> 


32 


(G 


1) 


2  _ 

Theoretical  estimates  indicate  G  ~  2  to  3 . 

From  a  survey  of  the  ultra-short  electron  beamlet 
simulation  results, ^  it  was  estimated  that  for  stable 
propagation 

g2  ,  J-EL— 

B  ttu  (1-r)  3 
0 

Since  the  physics  of  the  saturated  state  dees  not 
depend  on  beam  length  as  the  EM  pulse  is  many  synchrotron 
periods  in  length,  this  inequality  should  be  applicable  to 
finite  length  electron  beamlets,  and  stable  propagation  is 
obtained  by  choosing  S  =  l/v  large  enouah  to  satisfy 
Equation  (  2P)  . 


50. 


(  28  ; 
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However,  it  is  also  necessary  that  the  choice  of  £ 
to  secure  stable  propagation  be  in  a  parameter  regime  where 
the  linear  gain  is  finite. 

In  the  limit  of  a  long  electron  beamlet  with  no 
thermal  energy  spread,  the  linear  gain  tends  to  a  non-zero 
value  (Equation  27): 


Re  ( 1  -  r  +  6 )  K 


V-  3, 

4it  u  r  2 
o 


as  6  increases. 


1  >> 


ttu  62r5^  ru2 


3, 

Thus,  if  we  substitute  B2  =  nu0/2tT  2  in 
Equation  (28),  and  we  demand  Re  6  >  0 ,  we  obtain  the 
followina  inequalities  to  be  satisfied  in  order  to  have 
both  stable  propagation  and  finite  linear  gain: 


:<i-r)-uo 


21.5 


(29) 


T  (l-r)'tu 


2.26  r 


(30) 


These  ineaualities  impose  a  lower  limit  on  the  peak 
beam  current.  It  should  be  noted  that  in  deriving  these 
i  r.enua  1  i  t  i  es  ,  a  Lorer.tzian  beam  profile  was  considered 
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>  5.4,  this  lower 


(h(u)  =  l/~u„  (l+u'/ujl) .  For  r ^ u 
limit  is  determined  by  Equation  (30),  the  criterion  for 
finite  linear  gain.  This  implies  that  the  peak  beam  current 
I  must  exceed  a  value  given  by 


(1  -  r) 

k  r 
w  o 


(1  +  aw2)^ 


amps  (31] 


w 


where  a  =  eE  ''mk  c2,  the  wiggler  length  is  L  =  k  r2, 
and  rQ  is  the  beam  radius. 

The  corresponding  inequality  for  ultra-short  beamlets 
i  s 

— n —  >  _ L__  02) 

(1-r) f  0.06 

and  it  may  be  viewed  as  an  inequality  imposing  a  lower  limit 
on  the  beam  current  averaged  over  a  slippage  distance. 

Equation  (30)  imposes  a  current  limit  comparable  to 
Equation  (32)  for  7  ~  7.4. 

k. 

In  the  proposed  Los  Alamos  FEL  experiment  7  ~  6, 

while  for  the  nominal  parameters  discussed  in  Reference  2 
W 

.  ‘  ~  45.  Thus,  Equation  (31)  implies  a  required  peak  beam 
current  for  the  latter  case  larger  by  a  factor  of  ~  6  than 
that  discussed  in  Reference  (2). 

We  have  not  discussed  herein  the  question  of 
transition  from  linear  to  saturated  states  and  whether  any 
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difficulty  in  intermediate  states  arises  from  amplitude 
dependent  frequency  response,  i.e.,  chirping.  This  appears 
best  resolved  numerically  and  for  the  LASL  parameters  our 
simulations  confirm  behavior  in  linear  and  saturated  states 
with  a  smooth  transition.  Higher  efficiency  systems  (e.g., 

T  =  1000)  for  long  beamlets  require  very  lengthy  simulations 
which  remain  to  be  done,  although  for  ultra-short  beamlets 
(Reference  3)  no  new  difficulties  were  observed. 
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APPENDIX 


B 


SIMULATION  ALGORITHM  FOR  AN 
ARBITRARILY  TAPERED  WIGGLER  FEL 

A  relatively  simple  model  of  the  tapered  wiggler  FEL 
has  proven  adequate  for  the  analysis  and  simulations  involved 
in  determining  FEL  scaling  laws . ^  The  crucial  simplification 
involved  introduction  of  a  fictitious  uniform  static  electric 
field  into  the  model.  This  axial  field  served  to  maintain 
the  resonant  electrons  at  a  constant  y  as  energy  was 
extracted  by  the  ponderomoti ve  fields.  Thus  a  uniform 
wiggler  performs  as  a  linearly  tapered  wiggler. 

Recognizing  that  a  more  detailed  model  would  be 
necessary  to  optimize  and  evaluate  actual  tapered  wiggler 
designs,  the  following  model  was  derived.  It  has  the 
attractive  feature  of  retaining  the  basic  structure  of  the 
simpler  one  referenced  earlier,  thus  minimizing  the  recoding 
necessary  to  implement  it.  It  differs  from  the  previous 
model  also  in  that  the  optical  mode  is  assumed  to  be  a 
transverse  electric  rectangular  waveguide  mode  and  the 
wiggler  is  plane  polarized.  These  added  features  allow 
simulation  of  a  larger  class  of  physical  designs. 
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The  vector  potential  for  a  plane  polarized  wiggler  and 


transverse  electric  (TE)  optical  mode  were  taken  to  be 


A  =  x  A  (Z)  cos 


(K 


(  z )  dz 


x  Ag ( 2 , t )  cos  (kyy)  cos  (kgz  -  u)gt  +  £(z,t) 

(1) 

Here  kw,  Ag  and  C  are  the  wiggler  amplitude  and  wave 

number,  optical  amplitude  and  phase,  and  are  assumed  to  be 
slowly  varying  functions  of  their  arguments. 

Analysis  similar  to  that  of  Appendix  B  of  Reference  1 
yields  the  following  equations  for  the  particle  dynamics. 


d_ 

dz 


(Y  - 


V1 


It 

dz 


(Y  -  Yr) 


j  a  a 
s  w  s 

2y  c 
r 


sin  +  ;) 


(2) 

(3) 


dt 

dz 


(k  +  k  )  — 
w  s  ^ 

s 


where 


z 

..  =  Ik  dz  +  kz  -  ^  t 
J  w  s  s 

3 


eA 

w,  s 

a  =  - “ 

w,  s  me 


(4) 


(5) 

(6) 
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2c (k  +  5k  ) 
w  s 


(7) 


=  1  + 


w 


(8) 


5  k  =  k  -  — 
s  sc 


(9) 


Averages  have  been  taken  over  small-scale  oscillation  and 
has  been  assumed  that  y*  >>  aj  >>  as2  and  (Y-Yr)  <<  Y„ 
Also,  Y 2  >5>  1  and  3  <<  l  prior  to  entry  into  the 
wiggle..  All  of  the  particle  transverse  displacements  are 
assumed  negligible  compared  to  l/k^. 

The  transverse  velocity  may  be  computed  by  assuming 
conservation  of  transverse  canonical  momentum 


it 


10) 


Making  use  of  the  slow  variation  of  a  ^  and 
Maxwell's  eauations 


give 


j-A 


(  11  ) 


i  z 


c  ‘  k 


h)  %e‘;  ■ £ 


~e  a 

w 

k  cm', 
s  r 


•  x  . 


(  12  ) 
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where  H(z/V-t)  is  the  electron  pulse  number  density.  Note 

that  the  effective  photon  velocity  is  not  c  but  c2k  ju  , 

s  s 

which  is  the  group  velocity  of  the  TE  mode.  Coupling  to 
other  modes  has  been  neglected.  The  angle  brackets  represent 
averages  over  the  particle  phase  and  energy  distributions. 

The  filling  factor  f  is  defined  below  and  is  the  result  of 
averages  over  the  cos  k^y  dependence  of  the  TE  mode. 

Introducing  a  modified  definition  of  u  and  v,  see 
Equations  (18)  and  (19)  of  Appendix  C  of  the  1982  Annual 
Report , 
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Applying  this 
yields 
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to  the  equations  of  motion 
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beam 


(22) 
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M 
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is  the  total  number  of  electrons  per-  unit  area  in 
the  beam  pulse. 


The  functions  aw/a-  W'  -<w/  k -  w '  r^  {'r  an<^  y  which  are 
known  functions  of  z  must  be  transformed  into  functions  of 
u  +•  v  by  the  transformation 


u  +  v 


k  ‘ 

y  z 

k  k  L 

s  :w 

k  ; 
y 

k  k 
s  :w 


(23) 


Comparing  the  Equations  (15)  through  (17)  with 
Equations  (20)  through  (22)  in  Appendix  C  of  Reference  1,  it 
is  evident  that  the  algorithms  of  Equations  (26)  through  (52) 
of  Reference  1,  Appendix  C  may  be  utilized  with  the  introduc¬ 
tion  of  the  functions  aw/aoW.  kw/kjw,  - r/v  r  and  ky  as 
coe f  f ic i ents  . 

The  diagnostic  package  was  also  used  intact,  but  the 
interpretation  of  the  Fourier  transform  argument  within  the 
wiggler  must  include  the  revised  transformation  from  u,v  to 
z ,  t . 
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APPENDIX 


C 


SIDEBAND  INSTABILITIES 


The  linear  theory  of  sideband  instabilities,  involving 
the  coupling  of  small  amplitude  sideband  modes  to  the 
periodic  motion  of  electrons  trapped  in  the  ponderomoti ve 
potential  well  or  "bucket"  produced  by  the  combined  fields  of 
the  wiggler  and  a  large  amplitude  EM  pulse,  has  previously 
been  discussed  neglecting  perturbations  of  the  phase  of  the 
EM  pulse. *  More  recently,  the  analysis  has  been  extended  to 
include  both  amplitude  and  phase  perturbations  of  the  EM 
pulse . 

A  summary  of  the  principal  results  will  be  discussed  in 
this  appendix.  The  details  of  the  analysis  will  be  reported 
elsewhere . 


CONSTANT  PARAMETER  WIGGLER 

We  consider  stationary  states  in  which  the  ponderomotive 
potential  well  remains  constant  and  is  unaccelerated  through 
the  wiggler.  The  EM  pulse  is  represented  as  a  plane  wave 


N.  M.  Kroll,  P.  Morton  and  M.  N.  Rosenbluth,  "Free 
Electron  Lasers  with  Variable  Parameter  Wigglers,"  IEEE 
Journal  of  Quantum  Electronics,  Vol.  QE-17,  pp .  1436-1468 

(1981). 
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with  a  slowly  varying  amplitude  Ag  and  phase  Sg, 

Ag  exp  (-iwgt  +  ikgz  +  i£g).  The  one-dimensional  FEL 
equations,  written  in  terms  of  dimensionless  variables,  are 
described  in  Appendix  A. 

The  independent  coordinates  are  transformed  to  Z  and 
vs 


v 


c 


z 

c 


) 


where  L  is  the  wiggler  length  and  V  the  electron  resonant 
longitudinal  velocity.  The  front  of  the  wiggler  is  at  Z  =  0 
and  the  back  at  Z  =  1. 

The  linearized  FEL  equations  for  the  perturbed 
amplitude  Ag  =  A^g(Z)  e*<v  and  phase  =  ^(Z)  e  of 

sideband  modes  with  dimensionless  frequency  k  are  solved  to 
obtain  the  change  in  amplitude  and  phase  through  the  wiggler: 


a.  (Z=l)  -  a,  (0)  n 

- — ; — —  ■  r  Bi2 


(1) 


^(Z-l) 


CL(0) 


T-  a1C0)  Bn 


(2) 


where 


B12  *  2irJdJFo(J)  (t  *  l)  ‘  I1(K'V  <3> 
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J  =  —  j  E(k)  -  (1-k2)  K(k) 


*  I 


*/2 
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tt/2 


AL 


(1  -  k2  sin2  O  ' 


j  dc  (1  -  k2  sin2  O 

o 

exp  J  -  it  K(k 1 )  /K(k)  J 


=  (1-k2) 

2tt  f  dJ  F :  (J)  =  1 


(7) 


a0  is  the  dimensionless  amplitude  of  the  finite 
amplitude  pulse  and  is  independent  of  Z.  a^tZ)  is  the 
dimensionless  perturbed  amplitude.  nQ  is  proportional  to 

the  current  density  i/irr  2  ,  with  r  the  beam  radius. 

o  o 

f  is  the  filling  factor.  is  the  dimensionless 

frequency  of  oscillation  at  the  bottom  of  the  bucket. 
is  the  equilibrium  distribution  function  of  the  electrons 
trapped  in  the  bucket  and  is  a  function  of  the  electron 
action  variable  J  only.  Equation  (7)  determines  the 
parameter  k  as  a  function  of  J  or  vice  versa.  For 
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electrons  at  the  bottom  of  the  bucket  k  -»  0,  J  ♦  0,  while  at 

80o 

the  top  k  -*  1,  J  -*•  — — . 

It  is  assumed  that  n  /a  <<  1  so  that  the  linear  gain 

o  o 

per  pass  is  small. 

Let  a1(z»i)  *  e£a1(0)  *  (1  ♦  5  ♦  .  .  .)  a^O),  and 
ij^Z-1)  »  (1  +  5  +  .  .  •)  Cj(0)  where  5  is  the  linear  gain 
and  phase  shift  per  pass.  The  linear  gain  5  is  determined 
by  substitution  in  Equations  (1)  and  (2): 


[2.  |djro«)  (f  -  i)  -  ij 

[-2„{dJFo(J)  (f  -  l)  +  I2] 


(8) 


If  the  bounce  frequency  is  large,  >>  1,  gn  may  be 

approximated  by  gn  *  QQ/(ic-n^)#  and  1^  and  I2  reduces 

to 
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E  2*  n0  Jdj 


n  even 


SF  4  3  n  Q 

0  tt  n  q  0 

3J  K4  (1  -  qn)  2  (k  -  n.Q)  (10) 


In  the  limit  of  <  >>  Q  ,  l/(ic-nft)  may  be  expanded  as 

0 

a  power  series  in  </nP.#  and  by  doing  the  summation  in  ns 
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The  linear  gain  5  is  then  given  by: 
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The  magnitude  of  the  linear  gain  depends  on  the 
structure  of  the  electron  distribution  function  FQ(J). 

Two  specific  electron  distribution  functions  have  been 
investigated  and  the  linear  gain  evaluated  as  a  function  of 
the  sideband  frequency  K/^0' 


Case  I 

In  the  case  where  the  bucket  is  full: 
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J  <  J 


max 


(14) 
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The  integral 


•^nax 

2tt  J  dJ  Fo 

o 


1 

3 


and  the  summation  in  Equations  (9)  and  (10)  for  1^  and 

I_  ( Q  >>l  is  assumed  to  be  large)  may  be  converted  into 
2  o 

an  integral: 
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where  P 
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The  integrals  may  be  written  in  terms  of  4*-functions  or 
an  infinite  sum.  The  linear  gain  6  is  then  determined  by: 
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Re  _  is  tabulated  as  a  function  of  ■*—  in  Table  1 

n 

0  0 


Table  1 


< 

6a 

„  0 

Re  - 

0 

n 

0 

0.2 

0.059 

0.4 

0.196 

0.6 

0.332 

00 

• 

0 

0.412 

1.0 

0.428 

1.2 

0.393 

1.6 

0.265 

2.0 

0.147 

73 


No  comparison  with  simulations  has  been  made  because  of 
the  anticipated  difficulty  (imposed  by  running  time  limita¬ 
tions)  in  having  adequate  numbers  of  particles  very  close  to 
the  separatrix,  k2  -*•  1. 


Case  II 

In  the  case  where  F  decreases  linearly  in  k2  from 

o 

the  bottom  to  the  top  of  the  bucket 


64ft 


(1  -  k2).  1  >  k2  >  0 


(16) 


the  integrals  in  J  may  be  transformed  to  integrals  in 


/dJ  *  / 


dk 


4ft  K 


Thus  2tt  JdJFo(J)  |  -  l) 


0.6,  and 


the  other  integrals  may  be  evaluated  numerically.  In  the 

K2 

limit  <<  1,  Equation  (13)  yields  for  the  linear  gain: 


„  .  „  ''0  K 

0.62  —  — — 

a.  ft 

o  o 


(17) 


More  generally,  using  Equations  (5)  and  (6)  for  I 
and  I2,  the  linear  gain  6  may  be  calculated  from 
Equation  (8).  Ten  terms  in  each  sum  10  >  n  >  -10  are 
sufficient  to  obtain  convergence.  Re  6  is  tabulated  in 
Table  1  for  different  values  of  </ftQ.  The  values  of  the  FEL 
dimensionless  parameters  are  nQ  =  320  and  a^  =  640,  with 

ft  *  25.3  corresponding  to  four  bounce  periods  for  an 

o 

electron  at  the  bottom  of  the  bucket. 
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for 


w 


The  linear  gain  increases  linearly  with  </~ 

o 

<<  1  and  reaches  a  maximum  of  Re  5  ~  0.26  at 
</n0  =  0.75. 

The  growth  of  sideband  modes  has  been  simulated 
numerically  using  our  previously-described  I-D  partic'e- 
pushing  code. 

The  electrons  are  randomly  distributed  within  the 
trapped  region  of  phase  space,  weighted  by  a  linear  function 
of  k2  to  model  the  distribution  function  described  by 
Equation  (16).  The  procedure  for  determination  of  the 
sideband  eigenmodes  involves  three  simulation  runs,  one 
reference  run  with  only  the  main  optical  pulse  and  two  runs 
with  additional  linearly-independent  perturbations.  A 
response  matrix  analogous  to  Equations  (1)  and  (2)  is 
determined  by  Fourier  transforming  the  output  optical  pulse 
and  the  eigenvalue  of  this  matrix  yields  the  sideband  growth 
rates . 

The  linear  gain  evaluated  from  the  simulations  are  also 
tabulated  in  Table  2  for  different  values  of 

The  agreement  between  theory  and  simulations  is 
satisfactory . 

The  two  cases  discussed  in  this  appendix  are 
representative  of  the  classes  of  trapped  electron  distribu¬ 
tion  functions  which  could  be  found  in  realistic  FELs .  It 
may  be  noted  that  the  magnitude  of  the  linear  gain  of 
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Table 


0645  -i  0.0796  -i  0.0383  0.096  0.0511 


sideband  modes  as  well  as  the  structure  of  the  linear  gain 
dependence  on  frequency  for  these  two  cases  are  similar. 

This  analysis  can  readily  be  extended  to  discuss  the 
growth  of  sideband  modes  in  variable  parameter  wigglers  and 
is  currently  in  progress.  Some  preliminary  results  are 
presented  in  Figure  6  of  the  main  summary  of  this  report. 
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APPENDIX 


D 


PHASE  AREA  DISPLACEMENT 


I .  INTRODUCTION 


In  the  deceleration  of  relativistic  beam  electrons  by 
phase  area  displacement,^-  the  beam  electrons  are  injected 
into  a  variable  parameter  wiggler  in  which  the  resonant 
energy  of  the  wiggler  Yr  increases  from  the  front  of  the 
wiggler  to  the  back,  Yr(L)  >  Yr(0).  The  beam  energy  y 
is  such  as  to  produce  a  resonant  interaction  of  the  beam 
electrons  with  the  ponderomoti ve  potential  well  or  bucket 
(produced  by  the  combined  fields  of  the  wiggler  and 
electromagnetic  [EM]  pulse)  inside  the  wiggler,  y  =  Yr(z), 

L  >  z  >0,  usually  near  the  center  of  the  wiggler.  The 
interaction  may  be  viewed,  as  shown  in  Figure  1,  in  terms  of 
an  acceleration  of  the  bucket  through  the  phase  area  (y*^)  of 
the  beam  electrons.  The  result  is  a  downward  displacement  in 
energy  of  the  phase  area  occupied  by  the  beam  electrons  by  an 
amount  of  the  order  of  the  phase  area  of  the  bucket  divided 
by  2  it. 

If  the  bucket  acceleration  is  adiabatic  and  no 
electrons  are  trapped  in  the  bucket  during  the  interaction. 
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the  final  beam  energy  spread  is  nearly  equal  to  the  initial 
beam  energy  spread. 

This  method  of  beam  energy  extraction  has  two  features 
which  are  attractive  with  respect  to  operating  a  variable 
parameter  Free  Electron  Laser  (FEL)  oscillator  in  conjunction 
with  a  storage  ring:  (1)  The  energy  extracted  is  insensitive 
to  the  beam  energy  spread;  (2)  The  ratio  of  the  energy 
extracted  to  the  increase  in  the  root  mean  square  energy 
spread  can  be  made  to  be  small. 

In  practice,  the  bucket  acceleration  is  never 
completely  adiabatic  and  some  energy  spread  will  occur. 
Furthermore,  a  nonadiabatic  entry  and  exit  from  the  wiggler 
will  induce  additional  energy  spreads.  Thus,  such  a  device 
may  be  operated  in  steady  state  only  if  the  increase  in  the 
effective  phase  space  area  occupied  by  the  beam,  due  to  the 
energy  spreading  on  each  passage  through  the  wiggler,  can  be 
balanced  by  a  corresponding  decrease  due  to  incoherent 
synchrotron  radiation  in  the  storage  ring,  and  the  beam 
energy  boosted  to  compensate  for  the  losses  in  the  wiggler 
and  storage  ring.  A  convenient  measure  of  efficiency  is  the 
ratio  of  the  energy  radiated  in  the  wiggler  to  that  lost  in 
the  storage  ring  during  steady  state  operation.  It  is 
clearly  desirable  to  maximize  this  ratio. 

In  this  respect,  we  discuss  an  investigation  of  the 
efficiency  of  phase  area  displacement  wigglers  operated  in 
a  storage  ring.  The  startup  problem  is  first  addressed  in 


Section  I.  The  linear  gain  is  calculated  from  an  eigenmode 
analysis  of  the  linearized  FEL  equations.  This  is  followed 


in  Section  II  by  a  comparison  with  numerical  simulations 
of  the  theoretical  estimates  of  beam  energy  loss  and  beam 
energy  spread  due  to  interaction  of  electrons  with  a  finite 
amplitude  ponderomotive  potential  well.  In  Section  III,  the 
results  of  numerical  simulations  of  steady  state  operation  is 
described,  the  efficiency  is  determined,  and  the  scaling  of 
efficiency  with  the  wiggler  and  EM  pulse  parameters  is 
derived . 

A  full  self-consistent  simulation  would  be  very 
difficult  to  run  since  a  very  long  EM  pulse  is  required  for 
adiabaticity,  involving  a  correspondingly  long  electron 
micropulse  and  many  particles.  Our  results  must  thus  be 
interpreted  only  as  necessary  conditions  in  the  design  of 
phase  area  displacement  wigglers. 
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II 


LINEAR  EIGENMODE  ANALYSIS 


In  this  section  we  address  the  small  signal  problem; 
i.  e.,  growing  the  laser  pulse  from  a  very  small  amplitude. 

It  is  assumed  that  the  electron  beam  is  relativistic 
with  energy  parameter  y  >>'l  and  that  the  beam  self-fields 
are  negligibly  small  (Compton  regime).  The  periodic  magnetic 
field  wiggler  (with  wave  number  kw  and  length  L)  and  the 
EM  pulse  (with  frequency  wg  and  wave  number  kg)  are 
represented  by  circularly  polarized  vector  potentials  Aw 
and  Ag  with  negligible  spatial  variations  transverse  to  the 
direction  of  electron  beam  propagation.  The  EM  pulse  is 
approximated  by  a  plane  wave  Aa  exp  (-iw„t  +  ik^z  +  i£)  with 
slowly  varying  amplitude  Ag  and  phase  £ . 

The  electron  equations  of  motion  (in  the  resonant 

approximation)  in  terms  of  the  electron  energy  and  phase 

/\ 

angle  {y  ,4>)  and  the  equations  of  the  EM  pulse  amplitude  Ag 
and  phase  <;  (from  Maxwell's  equations)  are: 
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V  is  the  electron  longitudinal  velocity.  a  and  k  ,  and 

w  w 

hence  Yr ,  are  considered  to  be  constant  down  the  wiggler. 

The  actual  variable  parameter  wiggler  is  modeled  by  the 
inclusion  of  a  constant  electric  field,  that  is  7^0. 

Ayr  is  the  change  in  y  which  would  be  experienced  by  an 
electron  freely  accelerated  by  the  electric  field.  is  the 

nominal  electron  micropulse  beam  length.  h(u)  is  a  form 
factor  determined  by  the  beam  profile.  yq  is  the  beam 
radius  and  I  is  the  beam  current.  f  is  the  filling 
factor.  The  angular  brackets  imply  integration  over  the 
initial  energy  distribution  ( y  )  and  average  over  the  initial 
phase  (, _)  of  the  electrons: 


/ 


dy  F(ya) 
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F  (  y  )  is  the  initial  electron  distribution  function 

o 

normalized  so  that  I  d  i  F(y  )  =1. 

J  o  o 

In  the  (u,v)-space,  the  wiggler  (L  2  z  a  0)  lies 
between  the  lines  u  +  v  =  0  and  u  +  v  =  1.  The  electrons 
move  on  lines  of  constant  u  and  the  photons  of  the  EM 
pulse,  propagating  in  the  beam  direction,  move  on  lines  of 
constant  v.  The  beam  electrons  and  photons  interact  only 
when  their  trajectories  in  the  (u.v)-plane  intersect  within 
the  lines  u  +  v  =  0  and  u  +  v  =  1. 

The  electron  phase  space  (y,j)  trajectories  are 
determined  by  the  solutions  of  Equations  (1)  and  (2),  with 

initial  conditions  Y  =  Y  ,  'i>  =  V  at  v  =  -u.  The  phase  of 

oo 

the  electrons  entering  the  wiggler  is  uncorrelated  with  that 

of  the  EM  pulse.  Thus,  ;•  is  distributed  uniformly  between 

o 

0  and  2t. 

To  model  a  phase  area  displacement  wiggler,  a  constant 
decelerating  electric  field  (T  <  0)  is  applied  to  decelerate 
the  electrons.  At  the  front  of  the  wiggler,  the  electron 
energy  is  greater  than  the  resonant  energy,  y(0)  =  >  0. 

As  the  electrons  move  down  the  wiggler,  y  decreases,  passes 
through  zero  (at  which  point  the  electrons  are  in  resonance 
with  the  ponderomoti ve  potential  well),  and  is  negative  at 
the  back  of  the  wiggler  (the  beam  energy  is  less  than  the 
resonant  energy,  y(L)  =  Yf  <  o).  Thus,  the  deceleration  of 
the  electrons  is  used  to  model  the  increase  of  the  resonant 
energy  of  the  ponderomoti ve  potential  well  from  a  value  below 


the  beam  energy  at  the  front  of  the  wiggler  to  a  value  above 
the  beam  energy  at  the  back  of  the  wiggler. 

A  determination  of  the  linear  gain  from  an  eigenmode 
analysis  of  the  linearized  FEL  equations  has  previously  been 
discussed  in  the  1982  Annual  Report  for  FEL  oscillators  with 
constant  (7  =  0)  and  variable  parameter  wigglers  with  7  >  0. 
By  following  the  same  procedure,  the  linear  gain  for  FEL 
oscillators  with  a  phase  area  displacement  wiggler  (7  <  0) 
can  readily  be  derived. 

In  order  to  exploit  the  special  properties  of  a  phase 
area  displacement  wiggler,  it  is  desirable  to  have  an  EM 
pulse  with  frequency  such  that  the  electrons  (with  mean 
energy  y)  are  in  resonance  (with  the  ponderomotive  potential 
well)  near  the  middle  of  the  wiggler,  and  with  amplitude 
constant  over  many  slippage  distances  so  that  the  electron 
interaction  is  adiabatic.  In  order  to  grow  an  EM  pulse  with 
these  characteristics,  it  will  be  necessary  to  provide  for 
frequency  discrimination  (by  passing  the  pulse  through  a 
narrow  band  pass  filter)  to  select  the  mode  with  the  desired 
frequency,  and  to  use  long  electron  micropulses.  Thus,  it  is 
relevant  to  restrict  the  analysis  of  linear  gain  to  long, 
constant  amplitude  EM  pulses  and  correspondingly  long 
electron  micropulses.  In  this  limit,  the  emerging  EM  pulse 
is  a(L)  =  e',a(0),  where  the  linear  gain  and  phase  shift  per 
pass  5  (  |5|  <  1)  is  easily  derived  from  Equations  (1) 


through  (3)  to  be : 


Since  y  « implies  |e|  >>  1,  the  function  Z(V) 
may  be  approximated  by  its  asymptotic  expansions 


2(9)  =  2i  it  5  exp  {-&2) 


-  i  (1  +  1/2V  + 


. )  ,  Iro  &  <  0 


Thus,  if  Yfch  =  0, 


1  -  r  +  5  = 


2  r 


-  (1  +  i)  exp 


(— ) 

\  2 in  / 


/  i  ( I  r  I  -  y ) 2  \ 

+  (l-i)  exp  V — JTTl  / 


(6) 


The  linear  gain  oscillates  as  y  varies  and  will  be 
significantly  modified  by  a  finite  energy  spread. 

If  the  energy  spread  is  finite  with 


1  -  r  +  6 


i  n 

_ o _ 

2  |  r  !  e(e  -  i|r|*s) 
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and 


Re  ( 1  -  r  +  6 )  = 


noYth 


(27  -  I  r |) 


•\h  +  T’ 


(Yfc2h  +  (Y  -  H 


8n  y.  . 

— — —  (27/1  r  |  -  1)  ,  Y.  .  <  Y  a  I  r  1/2  . 

in3  th  (7) 


The  linear  gain  is  reduced  in  magnitude  and  is  positive 
only  when  y  >  \V\/2.  As  will  be  apparent  subsequently  when 
model  parameters  are  discussed,  the  energy  spreads  to  be 
expected  are  such  that  Equation  (7)  is  applicable  and  the 
predicted  linear  gain  is  too  small  to  be  useful.  Thus,  an 
alternate  startup  strategy  is  required;  e.g., 

(a)  Using  another  laser  for  startup. 

(b)  Varying  electron  energy  so  that  electrons  are 
initially  resonant  at  wiggler  entrance  or  exit. 
This  is  a  complex  scenario  and  it  is  not  clear 
that  it  is  feasible  when  the  transition  phase 
is  considered. 

In  the  remainder  of  this  report,  we  assume  that  this 
problem  has  somehow  been  solved  and  that  a  long  laser  pulse 
of  large  constant  amplitude  is  available. 
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III.  BEAM  ENERGY  LOSS  AND  ENERGY  SPREAD 

The  electron  equations  of  motion  [Equation  (1)  and 
Equation  (2)]  are  derivable  from  the  Hamiltonian 

/v  2 

H  ( Y » )  =  ~  -  r-|i  -  a  cos  (i}i  +  £) 

A  2 

=  Y  +  V(v)  (8 

Analytic  estimates  of  the  beam  energy  loss  and  beam 
energy  spread  can  be  obtained  in  the  limit  where  a  and  £ 
are  constant. 

In  this  limit,  5  is  a  constant  of  the  motion.  The 
electrons  see  a  potential  well  V(^)  »  -  Ttfr  -  a  cos  !j» 

(set  S  *  0),  which  is  plotted  in  Figure  2  for  T  <  0.  The 
stationary  points  of  V(v)  occur  at 

*  *  2Xt  ♦  ♦r,  2N~  ♦  vsign  (Vr)  -  $r 

where 

sin  *r  «  r/a 
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F(\|/) 


Electron  phase  trajectory  in  the 
pondTomotive  potential  of  a  phase 
area  displacement  wiaoler. 
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the  potential  well  depth  is; 


AV  =  2a  | cos  <l>r  +  sin  ^  sign  (<l>r))| 

corresponding  to  a  bucket  height  6W  of 

6^max  =  2a^  [cos  ^r  +  sin  'V  (^r  “  1  sign  (^r) 

The  difference  in  adjacent  maxima  is 

2  it  |  r  | 

2  it  a  sin  jij^  | 

A  typical  trajectory  in  the  ponderomoti ve  potential 
well  is  shown  in  Figure  2.  The  electron  phase  <P  increases 
from  an  initial  value  reaches  a  maximum  when  it  is 

"reflected,"  and  then  decreases  to  a  final  value  ij^.  Note 
that  if  the  amplitude  a  is  increasing,  the  electrons  may 
become  trapped.  Thus  the  wiggler  should  be  designed  to 
provide  a  very  slow  decrease  in  bucket  height.  Note  also 
that  the  amplitude  cannot  be  allowed  to  fluctuate  in  time, 
thus  requiring  a  very  strong  optical  filter. 

Thus,  by  integrating  Equations  (1)  and  (2)  through  the 
wiggler,  the  net  change  in  the  electron  energy  due  to 
interaction  with  the  ponderomoti ve  potential  well  is; 


Ay 

max 
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-Y  = 


y  -  rv 


v=l 

v=0 


TI 

■-/ 


a  sirO_d^_ 


^1 


^ 2(n-V(ip)j  J  yj2  f H  -  V ( -^ ) I 


J 


a  sin  d'l 


a  -  \[Y  a  sin  <l> 

“  s  r 


f 


di|j 


^  H  +  a  (\p  sin  i'  +  cos  i|/) 


V  H  +  a  (  iJj  sin  ip  +  cos  i  ) 


(9) 


where  Yf  is  determined  by  H  =*  -T'lp  -  a  cos  . 

The  integrals  in  Equation  (9)  have  been  evaluated  in 
Reference  1  and  the  following  estimates  obtained  of  the  beam 
energy  loss  and  root  mean  square  energy  spread: 


Ay  =  -  -  (a)  2  =  -  Ay. 

TT  L 


(10) 


wiggler 


sin  ^rl  Ayl 


(11) 
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A  further  energy  spread  can  occur  if  the  beam  entry 


into  and  beam  exit  out  of  the  wiggler  is  nonadiabatic .  This 
contribution  to  the  room  mean  square  energy  spread  is 
estimated  to  be 


/Ay  \ 

I  rmsl 


ends 


2a 

r 


(12) 


where  the  mean  initial  energy  7  ~  T/2.  However,  the 
magnitude  of  this  contribution  can  be  essentially  eliminated 
by  a  gradual  tapering  of  the  front  and  back  of  the  wiggler  so 
that  the  beam  enters  and  exits  the  wiggler  adiabatically . 

If  the  end  contributions  are  taken  to  be  statistically 
independent  of  the  wiggler  contributions,  the  total  root  mean 
square  energy  spread  is: 


Ay 


rms 


ay. 


(13) 


The  theoretical  estimates  of  A7l  and  Ay^, 

Equations  (10)  and  (13),  have  been  compared  with  the  results 
of  numerical  simulation  of  a  one-dimensional  FEL  with  a  phase 
area  displacement  wiggler. 

In  the  simulations,  a  monoenergetic  electron  beam  with 

A  A 

energy  y  =  y^,  distributed  uniformly  in  phase  angle,  is 
injected  into  the  wiggler.  The  transitions  into  and  out  of 
the  wiggler  were  step- functions ,  and  hence  nonadiabatic.  The 
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mean  energy  y ^  and  energy  spread  Ayrmg  on  exit  from  the 
wiggler  are  evaluated.  The  results  for  one  set  of  FEL 
parameters  (610  particles  are  used  in  the  simulations)  are 
shown  in  Table  1. 

It  may  be  noted  that  when  the  initial  beam  energy 

is  such  that  the  electrons  are  in  resonance  with 
the  ponderomotive  potential  well  inside  the  wiggler 

(T  -  2a^  >  >  2a  ;  ,  the  simulation  results  are  consistent 

with  the  theoretical  picture  of  an  average  energy  loss 
(A“  a  204)  and  an  average  energy  spread  (AYrms  a  40.0) 
independent  of  the  initial  energy,  although  the  simulated 
loss  was  about  20  percent  lower  and  the  spread  about 
25  percent  higher  than  predicted  by  the  crude  theory. 

No  significant  trapping  is  observed  except  when  the 
electrons  are  in  resonance  at  the  beginning  or  end  of  the 
wiggler.  When  trapping  occurs  at  the  beginning  of  the 
wiggler,  the  energy  spread  is  considerably  enhanced  as  the 
trapped  electrons  are  dragged  along  in  phase  space  by  the 
ponderomotive  potential  well. 


TABLE  1 


a 


104 


r  =  -1000 


Ay  =  Y*  -  Yf  ^rms 


150 

-212 

529 

-58 

436 

170 

212 

62.1 

200 

250 

214 

35.7 

207 

49.9 

300 

350 

198 

34.7 

400 

197 

34.6 

450 

202 

42.4 

500 

207 

30.7 

550 

205 

44.2 

600 

203 

38.1 

650 

203 

40.0 

212 

32.8 

700 

191 

39.1 

750 

203 

32.6 

800 

850 

204 

40.3 

203 

42.1 

900 

199 

58.7 

950 

1000 

196 

63.4 
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IV.  STORAGE  RING  OPERATION 


An  FEL  oscillator  with  a  phase  area  displacement 
wiggler  can  be  operated  in  steady  state  with  a  storage 
ring  if  the  increase  in  energy  spread  can  be  balanced  by  a 
decrease  due  to  incoherent  synchrotron  radiation  and  the  beam 
energy  boosted  by  a  radio-frequency  (RF)  cavity  to  compensate 
for  the  losses  in  the  wiggler  and  in  synchrotron  radiation. 

The  energy  loss  ^Syn  due  to  synchrotron  radiation 
is  proportional  to  Y2,  Ay  =  -S  y2.  In  terms  of  the 

syu  o 

dimensionless  variable  y: 


The  energy  change  per  pass  Ay  is  the  sum  of  three 


terms : 


= 


Ay  + 

syn 


1yrp  -  4yl 


*  —  Ay 


syn 


1  + 


Y  -  Y 

k  L 
w 


+  Ay 


RF 


'  Ay. 


(15) 


where  Ay^p  is  the  energy  increase  in  the  RF  cavity.  To 
attain  steady  state  operation,  Ay^  is  required  to 
compensate  for  the  total  mean  energy  losses,  ATSyn  +  Ay^. 

The  presence  of  a  storage  ring  can  therefore  be 
included  in  the  simulations  previously  described  by 
recirculating  the  electrons  through  the  wiggler  while 
modeling  the  effect  of  the  storage  ring  on  the  electron 
energy  with  the  following  equation: 


0  (n+1)  0  (n) 

Yi  =  o  y  +  b 


(16) 


where  y^n^  is  the  value  of  Y  on  exit  from  the  wiggler  on 
the  nth  pass,  and  -y^(n+l)  the  input  value  of  y  on  the 
(n+l)^  pass.  The  parameter  a  is  related  to  AySyn  by 


a 


Ay 


1  - 


syn 


k  L 
w 


(17) 
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The  parameter  b  is  determined  by  the  requirement  that  the 
mean  input  energy  of  the  beam  on  each  pass  is  maintained 
constant . 

The  simulations  are  carried  out  for  a  given  wiggler 
design  with  the  value  of  a  (and  hence  the  synchrotron  energy 
loss  AySyn)  adjusted  until  an  efficient  steady  state  opera¬ 
tion  is  obtained.  This  is  achieved  when  AySyn  is  just 
large  enough  to  limit  the  energy  spread  so  that  electron 
trapping  is  negligible.  Electron  trapping  is  undesirable, 
particularly  when  it  occurs  at  the  beginning  of  the  wiggler, 
since  it  results  in  considerably  enhanced  energy  spreading. 

The  evolution  of  the  electron  distribution  function  F 
per  pass  may  be  described  by  a  Fokker  Planck  equation,  and  in 
steady  state 


1 

2 


a2 

W 


<(Ay)2>  F  -  <  Ay  >  F  =  0 


(18) 


where  the  energy  spread  in  the  wiggler 

is  balanced  by  the  synchrotron  damping  in  the  storage  ring 


<  Ay  >  =  -  Ay 


(y  -  y) 

syn  k  L 
J  w 


If  the  dependence  of  (Ay rmg )2  on  y  is  neglected, 
the  solution  of  Equation  (18)  is: 
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F 


C  exp 
o 


k  L 
w 


(Y  -  7)  2 

(Ay  >2 

rms 


(19) 


where  CQ  is  a  normalization  constant. 

An  approximate  theoretical  criterion  for  negligible 
trapping  in  steady  state  may  therefore  be  obtained  by 
restricting  the  fraction  of  electrons  which  can  be  in 
resonance  with  the  ponderomoti ve  potential  well  at  the 
beginning  of  the  wiggler  to  less  than  10“^.  This  imposes  the 
following  lower  limit  on  A~s  n ; 


Ay 


syn 


k  L 
w 


log  10 


-4 


(Ay  ) 
rms 

L-  2 

(Y  -  2a V 


(20) 


where  ^  ~  |r|/2  so  that  electrons  with  the  mean  beam  energy 
are  in  resonance  near  the  center  of  the  wiggler. 

The  accuracy  of  this  criterion  has  been  estimated  by 
comparison  with  numerical  simulations  for  a  wiggler  design 
with  the  following  physical  parameters. 


B 


w 


k 


s 

L 


Peak  Circulating  Power 
Change  in  Resonant  Energy 


8  kG 
103 

4 

2^  x  10  cm 
2  x  10^  cm 
1.77  x  104 
4.9  GW/ cm2 
AYr  ■  28 
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The  corresponding  simulation  parameters  ares 

a  =  104 

| r |  =  103 

and  "y  was  set  equal  to 


Y  =  1 r 1/2. 

The  observed  minimum  value  of  Ay_„_/k,.L  for  which  no 

syn*  w 

significant  electron  trapping  occurred  after  200  steady  state 
passes  (using  123  particles  for  each  pass)  is 


Ay 


syn 


k  L 
w 


«  0.075 


(21) 


The  observed  energy  loss  in  the  wiggler  is 


&y 


laser 


»  205 


(22) 


and  hence 


Ay 


laser 


\ 


Ay 


/sii 


«  0.15 


(23) 


S^n  / simulation 

This  ratio  is  a  measure  of  the  FEL  steady  state  efficiency. 

The  theoretical  estimate  of  A"y  /k  L  obtained  from 

syn  w 

Equation  (20)  with  Ayrms  determined  by  Equation  (13)  and 


AyT  by  Equation  (10)  is 
Lt 
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Ay 

syn 

k  L 
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0.11 


(24) 


It  is  larger  than  the  observed  value  [Equation  (21)]  by  a 
factor  of  1.4.  The  difference  may  be  attributed  to  the 
dependence  of  Ayrms  on  y..  since  Aylager  »  AyL  is 
somewhat  larger  than  the  observed  value  [Equation  (22)],  the 
theoretical  estimate  of  FEL  efficiency 


Ay 


laser 


: 


Ay 


" 


«  0.13 


\  syn  /theory 

is  close  to  that  seen  in  the  simulations. 

The  scaling  of  efficiency  with  wiggler  and  EM  pulse 
parameters  may  be  obtained  from  Equation  (20)  by  substituting 
Equation  (11)  for  Ayrms  (the  energy  spreading  due  to 
nonadiabatic  entry  and  exit  is  assumed  negligible): 


„  k 

where  y  =  |r]/2  and  a  >  jr|  >  4a  . 

As  can  be  seen  from  the  preceding  discussion  of  the 
comparison  between  theory  and  simulation,  the  numerical 
coefficient  is  perhaps  uncertain  by  about  50  percent. 
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V.  DISCUSSION 


In  order  to  obtain  good  efficiency  during  steady  state 
operation,  it  is  essential  that  the  electron  interaction  with 
the  ponderomotive  potential  well  be  adiabatic  and  that 
electron  trapping  be  negligible. 

Furthermore,  a  long  constant  amplitude  EM  pulse  is 
necessary  to  produce  a  potential  well  which  remains  unchanged 
through  the  wiggler.  Variations  of  the  potential  well  can 
affect  not  only  adiabaticity,  but  lead  to  electron  trapping. 
As  a  general  rule,  the  electrons  should  "see"  negligible 
amplitude  fluctuations  in  a  traversal  time  across  a  "bucket." 
The  slippage  AZ  of  the  electrons  in  the  EM  pulse  during 
traversal  through  the  "bucket"  in  which  reflection  occurs 
(Figure  2)  is  estimated  to  be 


AZ 


>>  1 


excluding  those  exponentially  few  electrons  "caught"  at  the 
top  of  the  potential  well.  Thus,  an  approximate  criterion 
for  restricting  the  probability  of  trapping  to  less  than 
t  (due  to  increases  Aa  in  amplitude)  is  obtained  by 
demanding 

da 

- —  Az  ~  Aa  <  xa 
dz 
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that  is 


1  da 
a  dz 


t/  (2a) 


L  (?  - x) " in  (rfr) 


The  fraction  t  of  electrons  which  are  trapped  are 
presumably  accelerated  to  high  energies  and  must  be  replaced 
or  reprocessed. 

For  a  =  10^  and  r/a  =  10“^  and  r  =  10“^, 


1  da 
a  dz 


460L(^  -  l) 


which  implies  slow  amplitude  variations  over  many  slippage 
distances . 

For  the  numerical  example  discussed  earlier,  this 
implies  pulse  lengths  of  5  cm. 

When  the  above  conditions  are  satisfied,  the  FEL 
efficiency  in  steady  state  is  estimated  to  be 
[Equation  ( 25  )  ] : 


laser 

AT 


0.048  k 


(  P  (GW)  }  3'4  /  4a  2  \  (27) 

l*s!  t  l1'  in  I 

where  the  dimensionless  pulse  amplitude  j  a s  j  is  related  to 
the  power  P(GW)  in  gigawatts  of  the  EM  pulse  (assumed  to  be 
circularly  polarized)  by 
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P  (GW )  =  0 . 69  tt  k  2  r  2  a=  2 

S  p  5 


tt r  2  is  the  effective  area  of  the  pulse. 
P 


If  the  wiggler  length  L  is  taken  to  be  the  Rayleigh 


di f fraction 

length 

LR  =  ksrp  ' 

Equation 

(27)  reduces  to 

Ay, 

laser 

<  0.029 

(kwLR)5'  p 

/  a  2 

3/u  w 

f  fl  - 

AY 

syn 

3/, 

Y  z 

'  r 

\1  +  a  2 
w 

'  '  I r i '  (28 

The  extracted 

beam  energy 

a£(  ev) 

in  electron  volts 

per  pass  is 


A<?(ev)  = 


Ay 

laser 

0.511  x  10  Yr  ' 

w 


and  the  power  AP(GW)  transferred  per  pass  to  the  EM  pulse 
from  the  electron  micropulse  with  peak  current  I (amps)  in 
amperes  is 


AP (GW) 


-4  4a 

5.11  x  10  y  —r—r  I  (amps) 

L  IT  K 

W 


I (amps) 

(29) 
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In  steady  state,  the  power  lost  by  the  EM  pulse  at  the 
mirrors  equals  the  power  gained  from  the  electron  micropulse: 


AP(GW)  =  2(1  -  r)  P ( GW ) 


(30) 


and  hence  the  current  required  to  maintain  steady  state  for  a 
given  value  of  the  fractional  power  loss  2(1  -  r)  is: 


I (amps)  =  1 . 40  x  10' 


^7—  KV  )*  p,"(GKI 


=  1.4  x  10 


3  (1  -  r) 


(^)  pl/'<GW1 
w 


L  = 


k  r  - 
s  p 


(31) 


In  the  case  of  the  physical  parameters  described  in 
Section  IV  where  the  peak  circulating  power  is  P  =  4.9  GW, 
ks  a  2  7T  x  10^  cm-^,  y  =  10"*,  =  5.3,  if  the  power  loss  is 

1% ,  2(1  -  r)  =  0.01  and  irr  2  =  1  cm  ,  the  peak  electron 

kr 

micropulse  current  is  I  =  13.5  amperes. 

This  limit  on  the  peak  micropulse  current  is  not 
particularly  severe. 

A  more  severe  constraint  arises  if  the  operation  of  the 
FEL  oscillator  is  to  be  initiated  by  growing  the  desired  EM 
pulse  from  noise  levels.  The  peak  current  density  must  then 
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be  large  enough  to  produce  adequate  linear  gain  per  pass  to 
overcome  the  energy  losses  at  the  mirrors.  For  long  electron 
micropulses  with  finite  energy  spreads,  the  linear  gain  per 
pass  for  a  long  constant  amplitude  EM  pulse  is  determined  by 
Equation  ( 7 ) : 


k  2  L3  a  2  f 
—  J  w  w 

Re  S  a  3.76x10  — - - 5-  — - r—  I  (amps) 

Y  k  r  1  +  a  2 

r  s  0  \  w  I 


th 


(ft-) 


(32) 


If  the  energy  spread  is  taken  to  be  [Equation  (20)], 
the  maximum  allowable. 


Tth- 


|  Ay 


rms 


k  L 
w 


2Ay 


syn 


18.4 


the  linear  gain  for  the  example  previously  described,  where 
I  =  13.5  amps  and  the  beam  area  is  taken  to  be  itt^  =  1  cm2, 
i  s 

Re  5  =  1.1  x  10-3  -  lj  f  . 


This  value  of  the  linear  gain  is  somewhat  on  the  low  side  to 
overcome  any  reasonable  mirror  losses.  Higher  micropulse 
peak  currents  or  alternative  startup  strategies  would 
probably  be  required. 

Even  if  the  wiggler  and  beam  parameters  can  be  designed 
so  that  finite  linear  gain  is  obtained,  there  remains  the 
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issue  of  accessibility  to  the  desired  final  steady  state 
operation.  Further  investigations  are  needed  to  determine 
the  evolution  of  the  EM  pulse  from  noise  levels. 

In  summary,  high  FEL  efficiency  is  theoretically 
possible  for  FEL  oscillators  using  a  phase  area  displacement 
wiggler  in  conjunction  wit h  a  storage  ring.  Long  wigglers 
containing  many  wiggler  periods  and  EM  pulses  with  large 
circulating  peak  powers  are  effective  in  obtaining  high 
efficiencies.  Very  smooth  long  pulses  are  required  to  avoid 
trapping,  hence  good  frequency  discrimination.  The  linear 
gain  per  pass  tends  to  be  small  unless  appreciable  micropulse 
peak  current  densities  are  available.  The  scenario  of  the 
time  evolution  to  a  steady  state  is  complex  and  remains  to  be 
elucidated . 
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APPENDIX 


E 


TWO-DIMENSIONAL  EFFECTS  IN  FREE  ELECTRON  LASERS 


I.  INTRODUCTION 

The  motion  of  electrons  in  a  free  electron  laser  under 
the  combined  influence  of  the  wiggler  magnets  and 
electromagnetic  wave  may  be  calculated  from  the  particle 
Hamiltonian  as  discussed  below.  This  problem  has  hitherto 
been  considered  essentially  as  a  one-dimensional  problem  in 
which  the  electron's  motion  can  be  reduced  to  the  pendulum 
equation  in  the  "ponderomoti ve"  potential  formed  by  the 
wiggler  and  laser  magnetic  fields.1  In  particular,  high 
efficiency  "tapered  wiggler"  configurations  rely  on  electrons 
being  trapped  and  executing  "synchrotron"  oscillations  in 
this  potential. 

It  is  the  purpose  of  this  note  to  extend  that 
description  by  taking  account  of  the  fact  that  the  fields 
actually  depend  on  transverse  dimensions.  To  simplify  our 
discussion,  we  will  consider  transverse  dependence  on  only  a 
single  coordinate,  x.  The  most  obvious  consequence  of  the 
transverse  (focusing)  variation  of  the  wiggler  magnetic 
field  is  that  the  electron  undergoes  transverse  betatron 
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oscillations.  In  the  simplest  approximation,  the  betatron 
and  synchrotron  oscillations  are  decoupled,  since  one  is 
transverse  and  the  other  longitudinal  insofar  as  the  laser 
pulse  may  be  treated  as  a  plane  wave.  Herein  we  propose  to 
analyze  the  coupling  which  results  from  such  physical  effects 
as  curvature  of  the  optical  wave  front.  Radial  variation  of 
optical  amplitude  introduces  effects  similar  to  but  somewhat 
smaller  than  those  we  consider  here. 


Ill 


II.  DERIVATION  OF  DRIVEN  PENDULUM  EQUATION 


It  is  conventional  to  express  the  vector  potential  in 
dimensionless  units  a  =  eA/mc2.  We  denote  by  subscript  w 
quantities  relating  to  the  wiggler  magnet  and  by  subscript  s 
those  relating  to  the  laser  field. 


Hence  the  vector  potential  is 


Ayy  wi  th 


Ay  =  aw  cosh  kwx  cos  kwz 


-  a  (x,z,t)  cos  [k_(z  -  ct)  -  <j>(x,z,t)]  .  (1) 


Typically  aw  ~  eil)  and  ag  -  e^lO*3). 


The  transverse  dependence  of  the  laser  amplitude  and 
phase  are  presumed  given.  We  will  concentrate  our  attention 
on  phase  behavior  near  x  =  0  where  <t>  ~  ksx2/R  with  R 
the  effective  radius  of  curvature  of  the  wave  front.  As  a 
specific  example,  we  sometimes  consider  a  low  gain  oscillator 
for  which  the  laser  pulse  in  the  optical  cavity  will  be  well 
described  by  a  Gaussian  mode 

as  ~  (xo  +  )  *  exp  ('xVl  x0  *  2i  z/ks  I) 

'  s 
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the  width  of  the  beam  waist  and  z  the  distance 


with  x  „ 

o 

from  the  optical  waist.  For  such  a  pulse 


1 

R 


_1_  2/za 

2ZR  1  +  I2/zr)2 


(2) 


with  the  Rayleigh  length  ZR  =  ksxQ2  /2. 

We  may  use  the  canonical  momentum  pz  as  our 
Hamiltonian  with  (x,  px,  y,  t)  as  the  canonical  coordinates 
and  z  the  independent  variable. ^  Taking  y  >>  1,  kwx  <<  1, 
as  <<  1,  we  find 


Pz  =  Y  2y 


-  —  j  1  +  p  2  +  a2  [l  + 
2y  x  w 


k  2x2 

w 


cos  k  z 
w 


-  2a  a  cos  k  z  cos 
w  s  w 


k  (z  -  ct) 
s 


$ 


I 

I 


Averaging  over  the  fast  wiggler  oscillations  and  keeping  only 
the  term  in  the  cross  product  which  allows  for  near  electron 
resonance,  we  find  for  the  averaged  Hamiltonian 


Y  - 


|  1  +  p  2  +  (l  +  k  2  x2) 

2y  |  x  2  \  w  / 

(k  +k\z-kt-k  ~— 
\  w  si  s  s  R 


-  a  a  cos 
w  s 


(3) 
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Recalling  that  aw  <<  ag  we  find  for  the  transverse 
motion  the  decoupled  betatron  equation 


k  2  a  2 
w  w 


k6  x 


with  solution  x  =  x  cos  (k  Z  +6). 

3  3 

For  the  longitudinal  motion  we  have 


k  2  a  2  x; 
w  w 

2 


1  +  —  (1  +  k  2  x.2 
2  \  w  3 


k  a  a  r  , 

-  - •  - — -  sin  /k  +k  \z-  k  t  -k 
2y  w  s  I  s  s  R 


It  is  convenient  to  rewrite  Equation  (5)  in  terms  of 


the  optical  phase 


Ik  +  kiz-kct-k  — 
w  s  s  s  R 


l 


dz 


k  - 

w 


/  a 

k  (l  + 

(1  +  k  2  x  2)  ) 

s  \  2 

\  w  3  /  / 

2y ' 


k  x  2 

~R  k6Sin  (2V  +  26)  •  (8) 


Combining  Equations  (6)and  (8),  we  recover  the  familiar 
pendulum  equation  with  a  driving  term 


a  a  sin  . 
w  s 


2k  k,2  x  2 

s  3  B 

+  - - - -  cos  (2kg  z  +  26) 


or  on  introducing 


(2)  *  y‘ 


w 


1  +  —  (1  +  k  2  x  2) 
2  w  B 


(awas)  ; 


2k, 


z  *  £  z  ;  a  *  — —  ;  8  =  - 


,  k  x  2 

1  S_£_  ,2 


2  R 
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we  have  our  fundamental  equation 


j  2  .  , 

— pr  =  -  sin  v  +  3  sin  (az‘  +  2^)  .  (9) 

dz12 

The  phase  6  may  be  chosen  to  be  0  without  loss  of 
generality. 

Note  that  if  x.  is  comparable  to  the  laser  beam  width 

p 

xq  and  if  the  wiggler  length  is  of  order  a  Rayleigh  length, 
then  6  —©'(a2).  a  depends  on  many  parameters  but  typically 
a  ^  l,  depending  on  whether  the  peak  circulating  power  in  the 
laser  is  less  than  or  greater  than  1  GW.  These  limits  will 
be  refined  in  the  later  discussion,  but  we  anticipate  by 
noting  that  the  most  interesting  case  will  be  the  resonant 
case  a  ~  1,  and  that  such  a  power  level  is  likely  to  be 
realized  during  the  buildup  of  a  high  power  oscillator,  or 
possibly  in  the  operation  of  a  high  power  amplifier. 
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BEHAVIOR  OF  THE  BETATRON  SYNCHROTRON  RESONANCE 


We  now  turn  to  solution  of  the  basic  differential 
Equation  (9)  treating  ct  and  6  as  given  constants.  The 
solution  depends  on  whether  a  is  greater,  less  than,  or 
approximately  equal  to  1 .  We  bear  in  mind  that  the  primary 
interest  is  whether  particles  trapped  near  the  bottom  of  the 
ponderomoti ve  well  become  untrapped  due  to  the  driving  term. 


1 .  Case  I  -  ot  <<  1,  High  Optical  Power 

In  this  case  the  driving  term  changes  adiabatically 


3  sin ( ciz ’  +  2  6) 


Hence,  the  condition 


and  tp  «  sin'1 
tides 

a2  <1.  For  a  Gaussian  pulse  at  the  most  unfavorable 


that  particles  remain  trapped  is  simply  3  <  1,  i.e., 
ksx82 


value  Z  =  ZR,  we  see  from  Equation  (2)  that  this  condition 


R  ; 

x  £  , 

becomes  \  ot  <  1. 


Since  one  must  have  x  _  <  x  in 

B  o 


order  that  the  electron  remain  in  the  optical  field,  it 
follows  that  no  significant  detrapping  occurs  in  this  regime. 


2 .  Case  II  -  a  >>  1,  Low  Optical  Power 

Here,  since  we  are  interested  in  particles  with 
V  <<  1,  an  approximate  solution  is  ~  ^  J  sin  az'» 
leading  to  an  approximate  condition  for  avoiding  detrapping 
3/  :*2  <  tt/2  or  for  a  Gaussian  pulse  x;2/xq  <  2 tt  and  again  it 
would  seem  that  this  regime  is  not  dangerous.  The  form  of 
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the  solution  of  course  points  out  the  potential  dangers  of 
the  resonance  region  a  a  1 . 


3 .  Case  III  -  a  =  1 

Insofar  as  a  and  3  are  constants,  we  may  obtain  an 
analytic  approximate  solution  to  Equation  (9)  as  follows. 

The  Hamiltonian  for  Equation  (9)  is 


H  = 


-  cos  ^  -  hv  sin  otz 


If  we  introduce  angle  action  variables'^  for  the  unperturbed 
(3=0)  Hamiltonian,  we  have 


Here 


H  =  j  .':dJ  -  3 -i’  ( J,  $)  sin  az  ' 

J  =  JV  dv  =  ®  |E(k;)  -  (l-k:)  K(kz)j-, 


(10) 


Further 


and 


2K(kz) 


and  k 


2  _ 


H+l 


0  =  T 


r 

/ 


■y  =  4Kk/it 


J  k^  -  sin*1  i/2  / 


/  » f1?) 


7K  ( k 2 


=  4 


S  72^17  sin  (2n-^  V/cosh  (n-lz)"77PT 

n=l 

(ID 
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Here  E  and  K  are  the  usual  complete  elliptic  integrals, 
cn  is  the  Jacobian  elliptic  function,  and  K*  =  K  ( V 1 —X ^ }  . 
Electrons  trapped  at  the  bottom  of  the  well  are  represented 
by  k2  =  0.  Barely  trapped  electrons  have  k  =  1. 

Since  we  are  interested  in  the  case  of  near  resonance, 
we  substitute  Equation  (11)  into  Equation  (10),  keeping  only 
the  resonant  term  0  -  az '  .  Introducing  1  =  <J>  -  az'  by  a 
canonical  transformation,  we  obtain  the  new  autonomous 
Hami 1 tonian 


3 

2ka  -  a  -  E  -  (1  -  k2  )  K 


2  ;  co i 


cosh  {  tK ' /2K  f 


(12) 


which  is  now  a  constant  of  motion. 

The  bottom  of  the  potential  well  corresponds  to 
h'  =  0.  If  |cos  <t>  j  <  1  for  all  k2  between  zero  and  1, 
then  the  orbit  passing  through  the  bottom  of  the  well  is 
connected  to  the  top  and  most  particles  will  be  detrapped. 
Equation  (10)  is  easily  solved  for  cos<l>'  and  we  find 
graphically  that  such  detrapping  is  most  likely  to  occur  at 
a  =  0.865.  In  this  case,  3  <  0.1  is  found  to  lead  to 
detrapping.  Recall  that  we  have  previously  noted  that  for 
ot  «  1,  3  =  <  \  x  2/x  2.  Thus  total  detrapping  could  be 

DO 

avoided  for  x  2/x  2  <  0.4.  But  even  in  this  case 

B  o 

substantial  detrapping  could  occur  from  electrons  trapped 
away  from  the  bottom  of  the  well. 


A  more  important  limitation  on  predicted  detrapping 
results  from  the  finite  length  of  the  wiggler.  We  crudely 
estimate  the  detrapping  of  electrons  from  near  the  top  of  the 
well  where  7rK‘/2K  <  1.  From  Equation  (12), 


dJ 

dz 


dz'  3H 


=  2H  S 


sin  $ ' 


dz  3<J> '  's  cosh  ( 7t  K'/2K) 


(13) 


For  a  Gaussian  beam  using  Equations  (2)  and  (13)  and  putting 
Oi  k  1,  we  have 


Z/Z, 


•  2k~  _ — _ 

dz  »  xq2  1+  (Z/Zr)2 


(14) 


Those  electrons  with  J  nearly  equal  to  Jmax  -  3 / rr  may  now 
be  lifted  over  the  top  if  J  increases  sufficiently. 

Consider  a  wiggler  of  length  2L  with  optical  waist  at 
the  center.  Then  we  may  obtain  an  estimate  of  the  fraction 
of  trapped  phase  space,  AJ/Jmax,  which  moves  over  the  top  of 
the  well  by  integrating  Equation  (14)  from  Z  =  0  to  L.  We 
take  only  the  half  length  to  account  crudely  for  the  fact 
that  only  half  the  electrons  have  phases  leading  to 
dJ/dZ  >  0.  Then  the  possible  fractional  detrapping  is  given 
by : 


f 


AJ 


J 

max 


k 


S 


z 


R 


J.n 


(1+ 


(15) 
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Thus  it  would  appear  that  only  long  wigglers  with  large 
can  be  subject  to  serious  detrapping.  This  could  put 


an  important  emittance  constraint  on  long  wigglers.  However, 
we  note  that  another  such  constraint  is  already  conven¬ 
tionally  taken  into  account  in  wiggler  design.  This  arises 
from  the  fact  that  betatron  motion  at  fixed  y  leads  to  a 
change  in  electron  longitudinal  velocity  and  hence  does  not 
permit  electrons  to  be  trapped  in  the  ponderomoti ve  well. 

From  Equation  (8)  we  may  write 


L\ 

k 

/  k  a  \ 

s 

w  w 

\dz  ' 

P. 

\  2y  3  / 

v  'ef  f 

s 

< 


JT 


the  limit  coming  from  the  depth  of  the  potential  well  in 

Equation  (9).  This  means  that  the  maximum  allowable 

emittance,  ^  is  determined  from: 
ni  d  x 


k ,  Z 
3  R 


max 


=  2  si  2 


(16) 


Defining  the  maximum  allowable  emittance  from 
Equation  (16)  and  substituting  into  Equation  ((13),  we  have 
for  the  fraction  of  detrapped  phase  space: 


f  < 


4 


(17) 
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and  allowing 


For  the  usual  wiggler  L  =  ZR 
approximately  for  the  overestimates  we  have  made,  we  conclude 
f  ~  0.4  with  a  numerical  factor  uncertain  by  perhaps 

a  factor  of  2.  Pending  detailed  simulations  for  specific 
designs,  it  seems  clear  that  for  £/£max  <  %  the  detrapping 
is  not  significant. 
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IV.  SUMMARY  AND  CONCLUSIONS 


We  have  studied  the  effect  of  2-dimensional  motion  on 
electron  " synchrotron11  oscillations  in  the  ponderomoti ve 
potential  of  a  free  electron  laser.  We  find  that  the 
transverse  motion  of  the  electrons  is  essentially  an 
unconstrained  betatron  oscillation.  However,  if  the  optical 
wave  front  is  curved,  this  motion  couples  into  the 
synchrotron  oscillations.  We  find  this  coupling  to  be  weak 
except  near  resonance  where  2ka  *  ?.g  which  occurs  typically 
at  peak  circulating  optical  power  levels  near  1  GW.  Under 
resonant  conditions,  the  coupling  is  marginally  strong  enough 
to  lead  to  some  detrapping.  Thus  a  detailed  numerical 
simulation  of  specific  cases  may  be  required. 

However,  it  appears  that  in  an  amplifier  the  parameters 
may  change  rapidly  enough  that  the  resonance  is  passed 
through  without  significant  detrapping. ^  During  the  buildup 
phase  of  an  oscillator  passing  through  resonant  power  levels, 
we  have  roughly  estimated  peak  potential  detrapping  to  be  of 
order  of  40%  if  the  emittance  is  the  maximum  allowed  by  other 
considerations.  Hence,  if  a  reasonable  gain  margin  exists, 
this  detrapping  is  not  essential.  With  smaller  emittances, 
the  effect  is  of  course  smaller.  For  actual  cylindrical 
cases  where  many  orbits  do  not  pass  close  to  the  axis,  the 
effect  is  further  reduced. 
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It  should  be  noted  that  we  have  not  studied  the 
self-consistent  effect  of  betatron  motion  on  emission,  but 
only  the  effect  of  an  assumed  wave  shape  on  electron 
trapping.  With  these  caveats  we  conclude  that  2-dimensional 
effects  do  not  seriously  perturb  the  simple  1-dimensional 
picture  of  a  free  electron  laser  or  introduce  significant 
further  design  constraints. 

It  is  a  pleasure  to  acknowledge  discussions  with 
N.  Kroll ,  R.  Novick,  D.  Quimby,  V.  Wong  and,  particularly. 


D.  Prosnitz. 
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